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ABSTRACT 



Heat transfer effects of longitudinal vortices embedded within film 
cooled turbulent boundary layers on a flat plate were examined for 
freestream velocities of 10 m/s and 15 m/s, and for blowing ratios ranging 
from 0.47 to 1.26. Moderate strength vortices were employed having circu- 
lation to freestream velocity ratios of about 1.6 cm. Spatially resolved heat 
transfer measurements from a constant heat flux surface and mean tempe- 
rature distributions in spanwise planes show that local heat transfer is 
significantly affected by spanwise vortex position, and blowing ratio. Of 
particular significance are boundary layer and vortex structural changes 
which occur at high blowing ratios. 
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I. INTRODUCTION 



The increasing need for greater efficiency in gas turbine engines has re- 
sulted in higher turbine inlet temperatures. Consequently, combustors and 
turbine blading are subject to greater amounts of thermal stress, thermal fa- 
tigue, and creep. At present, gas turbines such as those associated with mi- 
litary applications, have inlet temperatures as high as 1800-2000 ^C, [Kef. 1). 
While the development of improved alloys with higher melting points is part 
of the solution, the development of efficient cooling systems is just as 
important [Ref. 2], In order to design an efficient cooling configuration, the 
heat transfer distributions for the gas turbine components are needed. 
Because of the complex geometries and flows involved near blades and end- 
walls, accurate convective heat transfer rates are difficult to obtain. [Ref. 31 
Ongoren, [Ref. 41, described the flow in a turbine cascade. As the inlet 
boundary layer approaches the blade, just in front of the blade, a horseshoe 
vortex forms. At the saddle point, it splits into a vortex on the suction side, 
and a vortex on the pressure side. The pressure side vortex becomes the 
passage vortex, moving from the leading edge of the blade towards the low 
pressure side of the adjacent blade. As the suction side vortex convects 
along the blade, it is eventually pushed away by the passage vortex from the 
adjacent blade. A smaller, corner vortex rotates in opposite direction to the 
passage vortex as was verified by Sieverding, [Ref. 51- 

The passage vortex is an example of a longitudinal vortex embedded in 
a film cooled turbulent boundary layer. In most cases, an embedded vortex 
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has a cross-siream length scale approximately equal to the boundary layer 
thickness. Therefore, the vortex is capable of strongly perturbing the 
boundary layer structure and modifying the heat transfer characteristics. In 
addition, longitudinal vortices usually maintain their coherence over a long 
streamwise distance, meaning that the heat transfer effects behind an effec- 
tive vortex generator are likely to be very persistent. (Ref. 6] 

Film cooling is used to provide a layer of cool fluid between a surface 
and high temperature free stream gases to which it is exposed. The film 
cooling not only acts as a thermal insulator but also as heat sink for the hot 
freestream gases. The overall effect of the film cooling is to reduce the 
temperature of the developing boundary layer, which in turn reduces the 
heat transfer to the surface. In the present study, the heat flux is calculated 
from : 

q = h(T w - T i*oo ) (eqn. 1 . 1 ) 

where h, includes the effect of the film cooling for constant wall heat flux, q ‘, 
and constant free stream temperature, T r co. In the present experiment 
variations of h correspond to variations of the wall temperature, T w . When 
T w = T r oo, equation 1.1 indicates that q" = 0, which may not be the case when 
film cooling is present. 

Goldstein and Chen, [Ref. 2], defined an adiabatic wall effectiveness for 
film cooling using: 



n = 




- T r 

- T r 



(eqn. 1.2) 
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The heat transfer is then calculated from : 



q = hf( Tv - Tav) (cqn. lo) 

Here, q "= 0 when T w = T aw with film cooling. Without film cooling T r oo = T av . 
and equations 1.1 and 1.3 are the same. 

Numerous studies have been conducted on the effects of film cooling on 
heat transfer, effects of secondary flows on heat transfer, and more recently, 
on effects of film cooling and secondary flows on heat transfer in a turbulent 
boundary layer. Summarizing the results of experimental research on the 
effects of secondary flows in turbine blades passages over the past decade, 
Sieverding (1984) concluded that it is absolutely essential to establish the 
significance of secondary flows at off design conditions, in particular in 
regard to leading edge vortices. 

Studying the effects of film cooling and secondary flows on heat trans- 
fer, Golstein and Chen (1983) performed an experimental study on the influ- 
ence of the end wall on film cooling of gas turbine blades using a single row 
of injection holes. Two years later, the same authors performed a similar 
study but this time employed two rows of injection holes for the film cooling, 
[Ref 7]. From these two studies, the authors concluded that, in the convex 
side of the blade there is a triangular region where coolant is swept away 
from the surface by the passage vortex, while the concave side was not sig- 
nificantly affected by secondary flows originating near the endwall. Other 
study on film cooling effects is given in [Ref. 81. 
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Eibeck and Eaton (1987), conducted a study of a longitudinal vortex em- 
bedded in a turbulent boundary layer over a constant heat flux plate. The 
authors found significant increases and decreases in local Stanton numbers, 
due to the thinning of the boundary layer on the downwash side of the 
vortex and thickening on the upwash side of the vortex. Spanwise heat 
transfer variations became larger as the circulation of embedded vortices 
increased. (Ref. 6] 

Investigating the interaction between a single weak streamwise vortex 
and a two-dimensional turbulent boundary layer, Russel, Pauley and Eaton 
(1987), observed a rapid growth of the vortex core, and a flattening of the 
core shape when the dimension of the core radius became comparable to the 
distance of the vortex center to the surface. Adverse pressure gradient 
caused an increase in the rate of core growth, and a stronger distortion of the 
core shape. These authors also provide an extensive review of work on 
vortices in boundary layers. (Ref 91 

In [Ref. 10], a study was conducted on heat transfer distributions and 
film cooling effectiveness on the endwall of an airfoil, using a full annular 
low aspect ratio vane cascade. The authors demonstrated the importance of 
the horseshoe vortices and secondary flows on the heat transfer and film 
cooling distributions. 

Recently, Joseph (1986), conducted a study of the effects of embedded 
vortices on heat transfer in film cooled turbulent boundary layers. The 
author showed that the effects of the vortex on heat transfer are significant 
and important : on the downwash side of the vortex, heat transfer is aug- 
mented, effects of the film cooling are negated and local hot spots will exist. 
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Near the upwash side of the vorte::, coolant is pushed to the side, appearing 
to augment the protection provided by the film cooling, (Ref. 3). Evans 
(1987), studied the fluid mechanics of vortices in boundary layers with film 
cooling, and showed that the vortices completely dominate the flow field, 
especially near the downwash side [Ref. 1 ]. 

The objective of this thesis is to increase the understanding of the effects 
on heat transfer, of a longitudinal vortex embedded in a film cooled turbu- 
lent boundary layer as affected by : 1 ) spanwise vortex position relative to 
the location of film cooling injection holes and, 2) blowing ratio. These 
effects are important regarding turbine blade and endwall heat transfer. 
The study was carried out under a series of steps. The first, was the design 
and construction of the heat transfer test surface, followed by qualification 
test to verify uniform wall heat flux, and an energy balance to identify and 
quantify the heat losses. After completion of the experimental apparatus, 
four types of heat transfer test were conducted : (1) heat transfer data with 
developing boundary layer only, (2) boundary layer and embedded vortex, 
(3) boundary layer with film cooling only, and (4) boundary layer with film 
cooling and embedded vortex. 
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II. EXPERIMENTAL APPARATUS 



A. WIND TUNNEL 

The wind tunnel pictured in Figure 3., built bv Aerolab, was designed to 
provide a flow field from the nozzle with uniform velocity and low turbu- 
lence intensity. It is designated the NPS Shear Layer Research Facility 
ISLRF). 

I . Description. 

The SLRF is a wind tunnel of the open circuit blower type with fan 
upstream and air entering the blower inlet from the surrounding room. The 
air speed through the test section can be adjusted from 3 to 40 m/s. The 
tunnel frame has leveling screws to adjust the center line of the tunnel to a 
horizontal position. The discharge of the fan slips into the inlet of the dif- 
fuser with a 1.6 mm clearance to isolate vibrations from the fan to the wind 
tunnel body. The diffuser section contains a filter pack and a nozzle leading 
to test section. The test section is a rectangular duct, 3.048 m long and 
0.6090 m wide. It is designed with numerous pressure tabs and four 38 x 
20.3 cm access ports along each of the side walls to provide easy access. The 
top wall is a continuous panel fabricated from 4.76 mm thick Lexan sheet, 
continuously sealed with neoprene along the edges. The ceiling height is 
adjustable to permit changes in the pressure gradient along the length of the 
test section. Additionally, the top wall contains numerous instrument ports 
for the measurements of various flow characteristics. The floor of the test 
section consists of three 0.6096 m and one 1.2 192 m long sections which are 
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removable and replaceable. These sections are all 0.6090 m wide and are 
sealed with "0" rings around the sides. Further discussions of the wind tun- 
nel are contained in [Ref. 1 1) and [Ref. 3: pp. 381. 

A schematic of the test section components with the coordinate 
system employed in the present study are shown in Figure l. An unheated 
starting length ot l 10 m exists upstream of the heated test surlace The in- 
lection nozzles are located 1.08 m downstream of the boundary layer trip 
and 0.02 m upstream of the test surface. The leading edges of the vortex 
generators are placed 0.479 m downstream of the boundary layer trip as 
shown in Figure 2. 

2. Qualification and Performance . 

Prior to its relocation, extensive qualification tests of the Shear 
Layer Research Facility were conducted by Ligrani, [Ref. 12). Results show 
that the variation of total pressure at the exit plane of the nozzle is less than 
0.4% at 26 m/s and 34 m/s. Mean velocity varies less than 0.7% for the 
same mean freestream speeds. From five-hole pressure probe measure- 
ments, the velocity angle deviation is nowhere greater than 0.6 degrees at 
the nozzle exit plane. 

Profile measurements of the mean velocity and longitudinal turbu- 
lence intensity in the turbulent boundary layer developing at 20 m/s indi- 
cate normal, spanwise uniform behavior. For this qualification test, and all 
results which follow, the boundary layer was tripped near the exit of the 
nozzle with a 1.5 mm high strip of tape. For the present study, total 
pressure measurements along the test section were made after relocation of 
the SLRF, in order to adjust the top wall for zero pressure gradient. The 
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maximum pressure difference along the test section was 0.007 in of H 2 O. 
Free stream turbulence intensity was measured to be 0.00085 i 8 5 one 
hundredths of one percent or 0.085%) at 20 m/s, increasing to 0.00095 at 30 
m/s. 

8. INJECTION SYSTEM. 

The injection system provides film coolant at temperatures above am- 
bient. The coolant is ejected from a single row of injection holes into the 
boundary layer developing along the bottom wall of the test section. The 
diameters of the injection holes were scaled relative to boundary layer 
thickness to be similar to a turbine blade, with 5/d ratio ranging from 0.37 to 
0.40. The free stream air is at ambient temperature, thus the direction of 
heat transfer is opposite that of a gas turbine. The temperature difference 
T w - T r oo range for this study has been kept less than 20 QC. to minimize the 
effects of variable fluid properties. 

The injection parameters M and 0 were scaled to resemble parameters 
near gas turbine blades where M ranges from 0.47 to 1.26 and 0 varied 
from 1.39 to 1.69. Since the average temperatures of the coolant, plate and 
freestream were approximately the same for different runs, T rc = 51 Q C, and 
T w = 40 Q C and T r oo ~ 20 Q C, variation of 0 was not a parameter considered in 
this study. Due to the reverse direction of heat transfer for the present ex- 
perimental apparatus T r c : T r <x> : T w ratio is 1.27 : 0.5- 0.55 :i.O as compared 
to 0.67- 0.83 : 1.5 : 1.0 for actual gas turbines. 
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1. Description 



Air for the inaction system originates in a two stage, 1 50 psig 
Ingersol-Rand air compressor, 10 HP, model # 71TD. From the compressor 
air flows through an adjustable pressure regulator, a cut-off valve, a 
reinforce flexible tubing (2.54 cm inside diameter), a moisture separator, a 
flow regulator, a Fisher and Portor rotometer (full scale 9.345E-3 m3/s, 19 8 
SCFM, model * 10A3565A), a diffuser, and finally into the injection heat 
exchanger and plenum chamber. The rotometer monitors the flow rate for 
film cooling. 

A photograph of the chamber, the row of film cooling holes and 
vortex generator is shown in Figure 3 b. The chamber is constructed of 1.27 
cm plexiglass, with outside dimensions of 0.305 x 0.508 x 0.457 m. The 
internal structure consists of three thin metal plates 0.381 x 0.508 m, 
starting 5 08 cm from the bottom and proceeding up at 5 08 cm intervals. 
Two silicon rubber heaters, 0.381 x 0.483 m, 120 volts, are separately placed 
over the bottom two plates. The heaters are controlled through a Powerstat 
variable autotransformer, type 136. The top surface contains 13 plexiglass 
injection tubes 8 cm long each, with an inside diameter of 0.952 cm i'3/8 ini, 
with a 1/d ratio of 8.42. The 13 injection nozzles are incline at angle of 30° , 
with a three diameter spanwise spacing between center lines. The middle 
tube is located on the center line of the test surface. 

Three pressure taps, positioned at the center of the front and two 
side faces of the injection plenum chamber, are used to measure P 0 c - Poo. 
Three 0.254 mm diameter copper-constantan wire thermocouples with 
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welded joints are placed at different locations inside the plenum chamber to 
measure T op . 

2. Qualification and Performance . 

Extensive qualification of the injection system was conducted by 
Joseph, [Ref. 3: pp. 23-26], Results show that the uniformity of the plenum 
chamber pressure, P oc . was satisfactory over the range of the injection con- 
ditions with differences of about 1% in the spanwise direction and a maxi- 
mum of 4% occurring for only one case at low flow rate of 0.327E-4 m?/s. 
This is equivalent of a blowing ratio M of 0.004 at Uoo = 10 m/5 

The plenum produces a reservoir of air at an elevated temperature 
and pressure, which is near stagnation conditions. The temperature at the 
exits of the injection tubes, T rc , is different from the temperature of the 
plenum of the injection chamber, T 0 p, due to conduction through the tubes 
surface to the surrounding air. It is more convenient to measure T op , 
whereas T rc is needed to calculate the injection parameters. Thus a 
relation between T rc and T op was needed and found from experiment to be 
Toe = 1-455 T op 0 86S . Because of the low velocities employed, (and negligible 
viscous dissipation), T oc ~ T rc within a fraction of a degree. 

In order to determine injection parameters, the following quantities 
must be measured, Poo, Too V c , P 0 c. T oc , and A, the area normal to the flow 
of the injection holes, designed to be 9.2633E-4 m2. The coolant velocity is 
then given by: 
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The static density is estimated using : 



Pco 

pc = w c 



(eqn. 2.2) 



The mass flux, m c , is the product of U c and p c . To calculate the isentropic 

mass flow. p c j and U c i are found using 

Poo 

pci = RT^ (eqn ' ' 



and 



(eqn. 2 A) 



The product of p c i and U c i for compressible flows may be given as : 




Discharge coefficients, Cd, are then estimated using : 

r PcUc 
d = (PcUc)i 



(eqn. 2.6) 



In order to check the injection system, discharge coefficients were 
measured as volumetric flow rate was changed. Results are given in Figure 
4 as a function of Reynolds number for 13 injection hole locations and 7 
injection hole locations. Discharge coefficients range between 0.5 and 0.730, 
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and increase with Reynolds number, for Re > 3 x HP, discharge coefficients 
collapse closely around a value of approximately 0.72 . Because these re- 
sults are consistent with those of (Ref. 8], satisfactory injection system 
performance is indicated. 

C. HEAT TRANSFER SURFACE. 

The heat transfer surface was designed and developed to provide a con- 
stant heat flux over its area. The average surface temperature may be 
adjusted and maintained from ambient up to about 60 Q C The plate was 
constructed so that its upward facing part is adjacent to the wind tunnel air 
stream, with minimal heat loss by conduction from the sides and beneath the 
test surface. The plate itself has been instrumented to measure tempera- 
tures with thermocouples placed just beneath the foil surface. A film of 
liquid crystals is sprayed over a layer of black paint painted on the foil. 

1. Description. 

A photograph of the heat transfer surface is shown in Figure 5 a. 
The design is based on ones used at the University of Minnesota [Ref. 1 3 i 
and [Ref. 14], and provides more uniform heat flux and better spatial resolu- 
tion of temperature than the surface used by Joseph, [Ref. 3: pp. 27-36], It 
consists of a thin stainless steel foil, AISI 302 full hard, 0.2032 mm x 1.3 m x 
0.467 m, painted flat black with 7 layers of liquid crystals. Attached to the 
under side of the surface are 126 copper-constantan thermocouples. 120 of 
these have flattened tips and were manufactured by Marchi associates (type 
MA 3%T-30-96-FEP-FJ), the remaining six have round junctions and were 
manufactured by Omega Engineering, Inc. Thermocouple lead wires are 
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located in grooves cut into a triple sheet ot' 0.254 mm (10 mil) thick double 
sided tape, manufactured by the 3M Company), as shown in Figure 5.b. The 
grooves are then filled with RTV. A thin foil constant heat flux heater. 1.0 
mm x 1.143 m x 0.457 m, 120 V/1500 W, manufactured by Marchi 
Associates, is attached to the tape with Electrobond epoxy. Beneath the 
heater is a 12.7 mm 11/2 in) thick lexan sheet, followed by a 25 4 mm of 
foam insulation, 82.55 mm thick styrofoam, three sheets ol 254 mm each 
lexan and one sheet of 9.53 mm thick balsa wood, as shown in Figure 7. 

Around the edges of the foil, grease was inserted between it and the 
plexiglass frame to fill any small gaps resulting from thermal expansion, tn 
addition, both the upstream edge and the trailing edge were taped to the 
bottom wall of the test section. Since additional vertical movement of the 
foil above the bottom wall of the test section occurs due to thermal expan- 
sion during heating, the level of the surface is adjustable and maintained 
level with the test surface by adjusting screws in the plexiglass frame sup- 
porting the heat transfer surface from below. During heat transfer tests, the 
top surface of the foil remained remarkably flat and smooth with minimal 
surface irregularities. The surface temperature is controlled by adjusting 
input voltage to the heater using a Standard Electrical Product Co. variac, 
type 3000B. With this type of heat transfer surface, good resolution of tem- 
perature was achieved without hot spots. 

Thermocouples are placed on the surface as shown in Figure 6. In 
each of the six rows, 21 thermocouples are located 1.27 cm apart from each 
other. The first row of thermocouples is located at 5 cm from the leading 
edge of the test surface, the second is 10 cm further downstream with 
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respect to the first, the third is 15 cm apart downstream from the second 
row, and the rest of the thermocouple rows are 20 cm apart from each other 
2. Qualification and Performance 

The present test surface was the second one constructed for this 
experiment. The first one could not be used since it had hot spots located 
under the thermocouple iead wires where the surrounding liner was not 
properly attach to the heater. In the second surface, this problem was 
remedied by increasing the thickness of the liner. 

To test the second heat transfer surface, a Huges Probeye Thermal 
Series 4000 Video System , consisting of a infrared and video camera with 
display screen, was used. Surface temperatures were measured as the test 
surface was heated outside the wind tunnel, since the system could not see 
through wind tunnel walls. With this system temperature variations as 
small as 1 2 C, can be measured. The surface was observed under three ope- 
rating conditions : natural convection with a surface temperature of ap- 
proximately 35 Q C, forced convection from the leading edge with a surface 
temperature of about 40 2 C, and force convection from the trailing edge, 
with a surface temperature of about 40 Q C. Results showed that most of the 
heat transfer surface temperatures were spanwise uniform within a fraction 
of a degree for all three tests. The only exceptions were several small cool 
spots with temperatures about 1 2 C lower than the rest of the plate, located 
between rows of thermocouple lead wires near the edges of the plate. 

A second test using the liquid crystals was performed to further 
qualify the test surface. The liquid crystals, manufactured by Davis Liquid 
Crystals, Inc., are rated from 30 to 35 2 C. Temperature variations as small as 
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1.2 can be measured with an uncertainty of t .2 y C The test was con- 
ducted with natural convection with a surface temperature of approximately 
32 Q C. Results were consistent with the infrared test with temperatures 
variations not greater than 1.2 2 C across most of the surface. As for the 
early test, small hot spots were evident near some thermocouple lead wire 
locations. 

3 . Energy B alance. 

An energy balance was performed to determine the heat loss by 
conduction from the heat transfer test surface used to obtain final results. 
During the energy balance, heat loss by radiation and convection were pre- 
vented since the metal foil surface, ordinarily exposed to convection in the 
wind tunnel, was covered with three layers of 25.4 mm thick foam insula- 
tion. For the energy balance, and for all wind tunnel testing, foam insulation 
was also placed around the sides of the test surface located below the wind 
tunnel convection surface. To estimate heat loss through the insulation 
placed on top of the foil surface, program ENERGB was employed with an 
algorithm involving the one dimensional, linear form of Fourier's conduction 
equation: 

AX 

q w = XA— (eqn. 2.7) 

For the insulation, K is 0.4 W/m C C, A is 0.4897 m 2 , aX is 0.0254 m, and AT 
is the temperature drop in the X direction in Q C. Heat conduction through the 
bottom and sides of the heat transfer device is then given by: 



q c - VI - q w 



(eqn. 2.8) 
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Here. VI is the power into the test plate, and q w is the conduction loss 
through the top insulation. Tests were made at five power leveis 11.02. 
13.0, 15.25, 16.52 and 18.2 Watts, chosen to give conduction losses at the 
same levels as experienced under normal operating conditions. Figure 8 



temperature. A second order polynomial was fitted to this data in order to 
predict conduction losses during heat transfer measurements: 



where Tdiff = T w - Tamb- This equation is valid over a range of T v - T a mb 
from 10 to 30 Q C. When exposed to convection in the wind tunnel, conduc- 
tion losses are only 1.5 to 2.5% of the total power, and thus, a 25% error in 
the estimate of conduction losses will cause less than a 0.5% error in the 
estimate of the heat transfer by convection. 

Radiation losses were estimated using two different approaches. 
For the first: 




q c = 0.683 + 0.954T diff - 0.0l6Tdiff2 



(eqn. 2.9) 




(eqn. 2.10) 



and 



n 

brad = Xdij 



(eqn. 2.1 1 ) 



j=l 
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[Rei. 151. For the second approach, 



n 




v Ji-Ji 

Zj(AiFij)- 1 



(eqn 2.12) 



[Rei. 16], The view factors, Fij, for the top and each of the side wails were 
estimated to be 0.54 and 0.23, respectively (Ref. 16], where 



From these two methods, the radiation losses for an average plate tempera- 
ture of 40 QC were estimated to be 55 Watts, approximately 8.5 % of the total 
power into the heat transfer surface. 

4. Contact Resistance Temperature Drop 

Local surface temperatures were measured using thermocouples 
placed in contact with the underside of the metal foil surface. The junctions 
of these thermocouples were held in place next to the foil with three layers 
of 0.254 mm thick double-sided lining tape and RTV silicon rubber epoxy. 
During the heat transfer tests, temperatures measured by the thermocouples 
were greater than those of the surface of the plate, due to thermal contact 
resistance between the thermocouples and the foil and conduction through 
the foil. The resulting temperature difference AT is given by: 



n 




(eqn. 2. * 2/ > 




( 1 AX\ 



(eqn 2.14) 
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where AX/kA accounts for conduction across the foii thickness and i/'h c A 
accounts for the thermal contact resistance, q is the heat flux through the 
foil. The contact resistance is highly dependent on the contact pressure as 
well as the area of the surfaces in contact. 

The value of contact resistance used in the present study was the same 
as used by joseph (Ref. 31. In his study liquid crystals were applied to the 
surface to measure its temperature distribution during convection tests, so 
that AT in eqn. 2.13 could be determined. Stanton numbers from baseline 
tests matched expected correlations after accounting for contact resistance. 
The empirical relationship for turbulent boundary layers at constant free 
stream velocity, along a flat plate with constant wall heat flux and unheated 
starting length of £ = 1.10 m. (Ref. 141 is given by 



St x p r 0A = o.030Re v "' 2 1- 



(£ 

w 



0.9 



-.1 1 1 



(eqn. 2.15. 



Data in Figures 10 and 13 show that this equation matches the data well for 
5 x 105< Re x < 2.0 x 10 6 . From these measurements, the two terms in 
brackets in equation 2. 13 were estimated to be 0.014 Q K/Watt. The same 
value was used for all thermocouples. However, on the present test plate, 
contact pressure and the area of the thermocouples in contact with the sur- 
face varied slightly from one thermocouple to another. Because contact 
resistance for each individual thermocouple may vary from this value, small 
deviations in the spanwise heat transfer coefficients and Stanton numbers 
result, which were independent of flow conditions. In the present study, the 
effects of these small variations were minimized by presenting results for 
local conditions in terms of Stanton number ratios . 
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D. TEMPERATURE MEASUREMENT. 



All thermocouples employed in the present study are type-T, copper- 
constantan thermocouples. 120 of these were manufactured by Marchi 
Associates, Inc and attached beneath of the heat transfer surface. These 
thermocouples, have a flattened junctions, so that they may be placed in 
good contact with the surface. One of these thermocouples was calibrated 
using a temperature regulated bath consisting of liquid nitrogen and electric 
heaters and a platinum resistance temperature reference (± 0.01 ^C). The 
calibration was performed over the temperature range of 0 - 45.32 Q C. A 
second order was used to convert microvolts to temperature : 

T - 0.018205 + 0.025846E - 0.000000581E 2 (eqn. 2.16) 



where: 



E = microvolts x 1 0 6 Volts (eqn. 2.17) 

This calibration was used for all Marchi Associates thermocouples attached 
to the plate, since all thermocouples indicated very similar behavior as 
temperatures were changed. 

Thermocouples manufactured by Omega Engineering Co. were used to 
measure six plate temperatures, as well as temneratures of the plenum air, 
freestream air, and the boundary layer as traverses were made. The calibra- 
tion of Joseph, (Ref. 3.: pp. 371, was first verified, and then used for all of 
these measurements except the traverses: 
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T = 26.573E - 1.937E2 + 0.998E3 - 0.26 1E^ 



(eqn. 2.18) 



Here, 



E = millivolts x 1000 



(eqn. 2.19) 



Temperatures surveys of T - Too, were performed using the automated 
traversing device. For this, two thermocouples manufactured by Omega 
Engineering were used, after calibration using the technique mentioned 
above. Simultaneous measurement of free stream temperature was made 
for every boundary layer location to minimize scatter due to wind tunnel 
temperature drift. For the thermocouple mounted in the probe that travels 
the flow field, the following equation was used to convert voltage to 
temperature: 



For the thermocouple used to measure the free stream, the equation 
employed was : 



As before, E = microvolts x 10 6 . The automated traversing mechanism has 
two degrees of freedom which allowed measurement of a plane the flow 
field. Each survey consisted of 800 probe locations, covering an area of 12 
cm x 22 cm. Both, spanwise and vertical traversing blocks are mounted on a 
20-thread per inch drive screw and two ground steel case hardened steel 
guide/support shafts. Each drive shaft is directly coupled to a SLO-SYN type 



T = 0.1836 + 0.025667E - 0.0000004882E2 



(eqn. 2.20) 



Too = 0.07832 +0.023 1054E - 0.0000006786E2 



(eqn. 2.21 ) 
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iV09 2 -FD3 1 0 steeping moior. Motors are controlled by a MIT AS two-axis 
Motion Controller. Both the motors and the controller are manufactured by 
The Superior Electric Co. The MITAS controller comes equipped with 2K 
bytes of memory and a MC68000, 16-bit microprocessor which allows the 
user to program the start, stop, duration, speed, acceleration and deceleration 
of the steeping motors. 

E. DATA ACQUISITION SYSTEM 

The data acquisition system was designed to rapidly measure thermo- 
couple voltages and convert them to temperature in degrees C. Using these 
temperatures along with user supplied information on ambient conditions, 
freestream conditions, power input, and flow rates into the injection 
chamber, the system calculates free stream velocity, density, local heat 
transfer coefficients, Stanton numbers, spanwise averaged Stanton numbers, 
and injection parameters such as blowing ratio and discharge coefficient. 

1. Hardware 

A Hewlett-Packard Series 300, Model 9836S computer, equipped 
with a MC68000, 8 MHz 16/32-bit processor, dual 5-1/4 inch floppy disk 
drives, and 1M bytes of memory RAM, was used in the study. The computer 
was used to process signals from the data acquisition system as well as the 
information supplied by the user in order to process, store, display and print 
results. An HP Think Jet printer was used to print data, and an HP 7470 two 
pen plotter was employed for graphics. Voltages from the thermocouples 
are read by an HP-3497A Data Acquisition/Control Unit with an HP-3498A 
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Extender. The unit communicates with the computer through a HP-829737 A 
Interface. 

2. Software. 

Programs STANFC1, STANFC2 and ACQTPRO were developed to 
process temperature and Stanton number data, program ENERGB was deve- 
loped to estimate the conduction losses, and programs PTSLC, PTSTAV, 
PLSTRVOR, PLSTRTIO, PLSTRFC, PLSTVV, SURFCONT and PLOTRUN were 
developed for plotting data. All these programs are listed in Appendix C. 
The data files created by the heat transfer programs are listed in Appendix 
D. 

Programs STANFC1 and STANFC2 are modified versions of STDAT1, 
which was developed by Ligrani, Ortiz and Joseph, (Ref. 3.: pp. 41], Program 
STAN'FCl prompts the user if film cooling is being used. If the answer is 
affirmative, the program prompts the user to enter the percentage of flow to 
the injection chamber from the rotometer and the pressure difference 
between the injection chamber and the static pressure in the wind tunnel. 
These inputs are transformed to SI units and stored in a data file named 
"FILDT". The program continues by prompting the user for the stagnation 
pressure of the free stream (in inches of H 2 O), the ambient pressure (in Hgj, 
the current (Amps) and voltage supplied to the heater. The program then 
reads the thermocouple voltages and converts them into temperatures. 
These are then stored in a data file named "TDAFC". After all temperatures 
have been calculated, the free stream density and velocity are calculated. 
Parameters pre iously calculated such as the free stream density, free 
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stream velocity, ambient pressure, and the power in are stored in a data file 
named "IDAFC". 

Program STANFC2, processes the data files IDAFC, IDAFC and FILDT. 
The program begins by prompting the user whether a vortex generator is 
being used or not, then it accounts for conduction, radiation losses and 
contact resistance. It continues by calculating local heat transfer coefficients 
and Stanton numbers, average spanwise Stanton numbers and Reynolds 
number based on downstream distance. These data are printed out and two 
data files, "STRFCV" and "STAY", are created. In "STRFC" the local heat 
transfer coefficients and Stanton numbers are stored, along with their 
spanwise position and downstream positions along the test plate. Spanwise 
averaged Stanton numbers with their corresponding down stream Reynolds 
numbers, are located in "STAV". The last section of the program includes a 
subroutine to calculate film cooling parameters, such as discharge coeffi- 
cients, density ratios, mass flux ratios, momentum flux ratios and blowing 
rate. 

Program ACQTPRO, developed by Ortiz, acquires temperatures from a 
thermocouple probe mounted on the automatic traversing device and from 
another thermocouple that senses the freestream temperature. This pro- 
gram was used to create the data to plot the temperature surveys of T - Tco. 
The program begins by prompting the user for the downstream distance 
from the boundary layer trip, the number of points in the spanwise direction 
where temperatures are going to be measured, the number of points in the 
vertical direction, the spanwise resolution in inches, the vertical resolution in 
inches and the initial position of the probe wtth respect to the plate, in both Z 



38 



and Y coordinates, in inches. A matrix of data points is then computed. 
Freestream temperature is calculated and ambient conditions are input. The 
program then enters a loop which samples each temperature 25 times per 
probe position and 5 times the free stream temperature. Upon acquiring 
these temperatures and obtaining their respective average, T - Too is calcu- 
lated for each probe location. The temperature difference aiong with their 
respective coordinates are stored in a matrix. At the end of the data 
collection run, these values are read into a data file named TPKO". on a 
floppy disk. These data are plotted using program "TPROPUN 
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Ill EXPERIMENTAL RESULTS 



A. BASE LINE MEASUREMENTS 

Heat transfer measurements were made at free stream velocities of 10 
m/s and 15 m/s, without vortex and without film cooling. These were done 
in order to determine how heat transfer from the plate compares with 
existing correlations. These measurements validated and qualified the per- 
formance of the heat transfer plate and measurements procedures used. 
Results are given in Figures 9-11 for 10 m/s, and in Figures 12-13 for 15 
m/s. 

Figure 9 shows the spanwise uniformity of the local heat transfer coeffi- 
cient of the test surface at 10 m/s. Except for row 1, the spanwise unifor- 
mity is very good with maximum variation of 10% (based on the average for 
a given row). These small variations are probably due to slight differences 
in the spatial uniformity of the heat transfer test surface, especially diffe- 
rences in the conduction contact resistance between different thermocouples 
and the metal foil, which comprises the top of the heat transfer plate. Larger 
spanwise variations for row l are believed to be due to multi-dimensional 
heat transfer by conduction through the leading edge and through the front 
corners of the test plate in addition to the contact resistance. These varia- 
tions are not as large as those observed by Joseph [Ref. 3: pp. 45, 77-78). 

The spanwise-averaged Stanton Numbers for 10 m/s are shown in Fi- 
gure 10. Six data points are shown, where each corresponds to a different 
thermocouple row. The data shows excellent agreement with the empirical 
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equation for turbulent boundary layers on a flat plate at constant free 
stream velocity, with constant heat flux and unheated starting length of x = 
1.10 m [Ref. 151- 



A survey of T - Too in degrees Centigrade was obtained at X = 1.48 m, 
using a calibrated copper-constantan thermocouple, positioned using the 
automated traversing device. Results are shown in Figure 11, for Uoo-10 
m/s with no vortex and no film cooling. In this figure, results for an area 
normal to the flow of 1 1 cm x 20 cm. are given, as measured using 800 
probe positions. The survey evidences good spanwise uniformity of the 
temperature boundary layer. 

Figure 12 and 13 show the spanwise variations of the heat transfer 
coefficient and averaged Stanton numbers as function of Reynolds number 
for 15 m/s free stream velocity. Results are consistent with those obtained 
for 10 m/s. In Figure 12, the data show good spanwise uniformity, except 
for the first row, where variations are slightly larger than for the 10 m/s 
case. In Figure 13, the data again show excellent agreement with equation 
(3.1). 

B. BOUNDARY LAYER WITH SINGLE VORTEX 

Heat transfer measurements were made at free stream velocities of 10 
m/s and 15 m/s with an embedded longitudinal vortex and no film cooling. 
These were done to further qualify the heat transfer surface for spatially re- 
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solved heat transfer measurements. To produce the vortex, vortex generator 
*2, (see Joseph (Ref. 3: pp 761), was positioned 0.479 m downstream from the 
boundary layer trip, and 4.79 cm in the positive Z-direction, as shown in 
Figure 2. 

The spanwise variation of local St/Sto for 10 m/s and 13 m/s are pre- 
sented in Figures 14 and 15, respectively. Here, Sto is the local Stanton num- 
ber without an embedded vortex. In both of these, the effect of the vortex 
on wall heat transfer is evident. The Stanton number ratios approach 1 
away from the vortex. As Z decreases from * 1 5cm, the ratios increase until 
maximum values of 1.2 and 1.25 are reached for 10 and 15 m/s, respec- 
tively. These maxima correspond to locations on downwash sides of vor- 
tices, where secondary flows result in boundary layers which are locally 
thinner than at other locations. For -5.0<Z<.0.0 cm a large gradient of heat 
transfer exists at each downstream location. The location of this gradient 
moves to smaller Z as the vortex develops downstream. For smaller Z, 
Stanton number ratios for 10 and 15 m/s, reach minima of 0.92 and 0.90 at 
locations which correspond to the upwash sides of vortices. Here, the 
boundary layer is locally thicker than at other locations. 

The results in Figures 14 and 15 are in excellent quantitative agree- 
ment with those of Joseph [Ref. 3: pp. 87-881, for the same experimental 
conditions, (The sign of the Z coordinate in this experiment is opposite to 
the one employed by Joseph). The results in 14 and 15 also show qualita- 
tive agreement with the measurements of Eibeck and Eaton (Ref. 6], where 
small quantitative differences are a result of different vortex generator 
geometries. 
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C. BOUNDARY LAYER WITH FILM COOLING 



In order to further qualify the heat transfer surface and the film cooling 
injection system, measurements were made with film cooling without em- 
bedded vortex. Results were obtained at 10 m/s , at 0 = 1.614, 1.699 and 
1.50 with blowing ratios of 0.68, 0.98 and 1.26 respectively, and at 15 m/s. 
at 0 = 1.395 and 1.586 with blowing ratios of 0.47 and 0.86 respectively 
Results are presented in Figures 16-30. 

Figures 16-18 show the spanwise variations of St/S to at 10 m/s, tor 
blowing ratios of 0.68, 0.98 and 1.26, respectively. The first data set was 
obtained using 13 injection holes, which produce a film wide enough (25.4 
cm) to cover the entire span of the measuring heat transfer surface. Because 
of the limitations on the secondary air supply, 9 and 7 injection holes were 
required to produce blowing ratios of 0.98 and 1.26, respectively. In these 
cases the coolant was injected such that it covered a portion of the center 
span of the heat transfer test surface. Consequently, Stanton number ratios 
are higher near the edges of the plate, as shown in Figures 17 and 18 The 
latter Figure for only 7 injection holes, shows this effect to be particularly 
significant, however, in spite of this a large spanwise uniform area exists in 
the center portion of the plate where adequate spanwise averages may be 
obtained. Apart from this, the spanwise uniformity of Stanton number ratios 
in Figures 16-18 is excellent with a maximum deviation of about 10 percent. 

Figures 19 and 20 show the spanwise variations of the Stanton number 
ratios at 15 m/s. Results for blowing ratios of 0.47 and 0.86 were obtained 
using 13 and 7 injection holes respectively. Spanwise variations in Figures 
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19 and 20 are larger than in Figures 16-18, though qualitative trends are 
very similar. 

A summary of spanwise averaged St/Sto as dependent upon Reynolds 
number is given in Figure 21. Graphs showing individual curves are pre- 
sented in Figures 22- 26, where averaged St/Sto are given for blowing 
ratios of 0.68, 0.98 and 1.26 at 10 m/s, and 0.47 and 0.86 at 15 m/s, respec- 
tively. Reterring to Figure 2 1, the lowest St/Sto are observed for M = 0 68 at 
10 m/s and for M = 0.47 at 15 m/s. For each individual data set, the lowest 
ratio is present for the two rows of thermocouples nearest the film cooling 
holes, where X/d is 7.35 and 17.48. Here, X is distance from the down 
stream edge of the injection holes, and d the injection hole diameter. The 
data in figures 21-26 are plotted as function of the blowing ratio in figures 
27 and 28. Results in 21-28 are qualitatively consistent with the fully tur- 
bulent measurements of Goldstein and Yoshida, [Ref. 151, also for a single 
row of injection holes. 

Two surveys of T - Too in degrees C were obtained in the film cooled 
turbulent boundary layer at X = 1.48 m, X/d = 41.94. Both surveys were 
made at a free stream velocity of 10 m/s, with M = 0.98 and a coolant injec- 
tion temperature of 5 1 Q C. In one case, the plate was heated to 40 C C (Figure 
30.), and in the other, the plate was maintained at the free stream tempera- 
ture (Figure 29 ). The latter experimental arrangement with unheated plate 
allowed the determination of the presence and distribution of heated film 
coolant as indicated by higher temperatures. In Figure 30 with a heated 
plate, the locations of coolant jets are evident as local hot spots near the wall. 
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Comparing these results to ones in Figure 1 1. without film cooling shows that 
the thermal boundary layer is about twice as thick. 

D. BOUNDARY LAYER WITH SINGLE VORTEX AND FILM COOLING 

In this part of the study the effect of an embedded vortex on a film 
cooled turbulent boundary layer was studied. This section of the experiment 
was conducted in four different parts. For all the parts of the experiment 
the same vortex generator was employed. The first consisted in repeating 
Joseph data, [Ref. 3. : pp. 97-104): Uoo= 10 m/s, vortex generator # 2 at Z = 
4.79 cm and blowing ratio, M = 0.68, in order to provide additional check on 
measurements, equipment and procedures. In the second step, a blowing 
ratio of M = 0.98 was employed, at Uoo = 10 m/s. The position of the vortex 
was changed to three different locations with respect to test section center 
line. Vortex position A corresponds to vortex generator at Z = 3.52 cm, vor- 
tex position B to Z = 4.79 cm, and vortex position C to Z = 6.06 cm. In the 
third step, a freestream velocity of 15 m/s, with M - 0.47 were employed, 
in order to study the effect of blowing ratio. In the fourth part, a blowing 
ratio of M = 1.26, at Uoo = 10 m/s was employed, in order to study the effect 
of high blowing ratio. Here, information was also obtained on changes in 
boundary layer structure which occurred for higher blowing ratios. 

1. Freestream 10 m/s. Bloving ratio. M = 0.68. Reneat of 
Joseph data. 

These heat transfer measurements were obtained to compare with 
those of Joseph, (Ref. 31. He first showed that embedded vortices cause sig- 
nificant changes in heat transfer in film cooled turbulent boundary layers. 
This finding is also evident from the present study. To obtain the present 
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data set, vortex generator # 2 was placed at position B, and all 15 injection 
holes were used tor film cooling. Results are presented in Figures 3la.-31.f, 
and in Figure 32. Overall effects are now discussed. The embedded vortex 
produces a thicker boundary layer in the upwash side thus augmenting the 
film cooling protection. Near the downwash side, the boundary [aver is lo- 
cally thinner, and the protection provided by the film cooling is minimized. 
The undercooled region produced on the upwash side vortex is very per- 
sistent not only in the down stream direction, but also in the spanwise 
direction. Here, Stanton numbers are .5 to 5.5 percent lower than those in 
boundary layers with film cooling only. 

Figures 31.a-f„ show the spanwise variation of the Stanton number 
ratios, St/Sto with film cooling and embedded vortex. Here, St is the local 
Stanton number and Sto is the Stanton number for the same location and 
free stream velocity but without film cooling and without embedded vortex. 
A St/Sto = 1 thus indicates an undisturbed thermal boundary layer. As 
shown in the set of Figures 3 1 .a - 3 If, the lowest St/Sto is 0.55 at X = 1.15 m, 
X/d = 7.35, and the maximum is 1.025 at X = 2.00 m, X/d = 96.50. If these 
results are compared with those of Joseph (Ref. 3. : pp. 97-1031, a general 
qualitative agreement can be observed. Differences are due to im- 
provements in the present heat transfer test plate, especially better 
spanwise temperature resolution and more uniform heat flux. 

Figure 32. shows surface contours of St/Stf. Here, St is the local 
Stanton number, and Stf is the local Stanton number at the same location 
with film cooling but without embedded vortex. A very steep heat transfer 
gradient is present which is near the same location as the axis of the vortex. 
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At the upstream end of the plate this gradient is located near Z = - 2 cm. A 
hot spot is present at larger Z, and a region of high heat transfer is present at 
smaller Z on the upwash side of the vortex. These results are in qualitative 
agreement with those of Joseph [Ref. 3.: pp. 104], except the most significant 
hot spots are located further downstream in the present study. 

2 . Fr ee stream 10 m/s. Blow ing ra tio . M ^ 0 .98 

In order to achieve the blowing ratio used for this data. 9 injection 
holes were employed for film cooling. Here, injection covers the midspan of 
the test plate from -12 cm < Z < 12 cm. Vortex generator *2. was placed 60 
cm upstream of the row of injection holes, at vortex position A, B and C to 
investigate the effect of the vortex position with respect to film cooling 
injection holes. Results are presented in Figures 3.3 - 42. 
a. Data overview. 

Spanwise variations of St/Sto for the vortex position B are 
shown in Figures 33a-33f. These data show the same overall qualitative 
trends as measurements at 10 m/s and M = 0.68 (Figure 31). With embed- 
ded vortex, the normalized Stanton numbers in 33a-33f are increased on the 
downwash side and decreased on the upwash side relative to the boundary 
layer with film cooling only. The only significant quantitative difference 
between these results and those in 31 are higher St/Sto at X/d = 17.48 and 
at X/d = 33. 39 resulting from a less effective film cooling at M = 0.98. 
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b. Effect of downstream development. 

In Figures 33a - 33f, a double peak of St/Sto is present at X/d 
= 17.48, and at X/d = 33.59. As the boundary layer develops further down- 
stream, one peak becomes smaller (Z ~ -2.0 cm), and one increases in mag- 






1.0 cm), until at X/cl = 54.59 the St/Sto distribution shows ora. 



iarge peak. This St/Sto peak increases in magnitude to become equal to 1.05. 
or about 34 percent greater than the St/Sto with film cooling only. St/Su' 
data in Figure 34. indicate that this hot spot cov rs an area along the center 
and spanwise downstream half of the plate. Such behavior evidences highly 
three-dimensional interactions within the boundary layer. 

Figure 34 shows the St/Stf distribution for vortex position B 
and M = 0.98. A steep heat transfer gradient is evident along the the lengtn 
of the test surface, which corresponds to the path of the vortex center. This 
path is at an angle to the streamwise direction, moving to smaller Z as the 
vortex convects downstream. Compared to results for M = 0.68 in Figure 32, 
the heat transfer gradient is steeper, and a larger region where St/Stf is less 
than 1 is present over the top third of the test plate. This follows since a 
larger amount of coolant is swept by the vortex at higher blowing ratio. 

For Z values smaller than those along the location of the vortex 
center, Figure 33 shows values of Stanton ratios which are about 5 per cent 
lower than those with film cooling only. The same trend is evident in Figure 
34. where almost 1/3 of the plate shows 0.90 < St/Stf < 0.98. These locations 
are on the upwash side of the vortex, where film coolant seems to be pushed 
and spread over a very iarge portion of the plate. Similar phenomena were 
observed by Goldstein and Chen. [Refs. 2 and 71, in a study of film cooling of 
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a turbine blade with injection through one and two rows of holes in the 
near-endwall region. 

Figures 35a- 35d show the T - Tco temperature field as it de- 
velops downstream at X/d = 41.9, 82.9, 109.2 and 147.0. For these tests, 
film cooling jets were heated to 51 2 C without providing anv heat to the test 
plate. Thus, the temperature field shows how fluid from the film cooling 
holes is convected and distorted by the vortex, where higher temperatures 
indicate greater amounts of coolant. The most dramatic effect evident from 
Figures 35a-d occurs on the upwash side of the vortex, where coolant is 
lifted away from the thermal boundary layer which is ordinarily about 4 cm 
thick with film cooling. At X/d = 41.9 fluid affected by the film cooling is 6.5 
cm from the surface, and at X/d = 147.0, the temperature field is more than 
10 cm from the surface. Also, evident are cooler temperatures with down- 
stream development resulting as coolant is convected and diffused. The 
areas indicated by low St/Sto and St/Stf values correspond to 0. < Z < -5. cm, 
where coolant seems lifted from the surface to accumulate in one small area. 
Figures 35a-d clearly show how the coolant is pushed and spread over an 
area by the vortex, and how this area increases in size, as the flow convects 
downstream. 

On the downwash side of the vortex, secondary flows cause 
freestream fluid to be located very near the wall. Here, the thermal bound- 
ary layer is greatly thinned and very little effect of the film cooling remains. 
From Figure 35c for X/d = 82.94, the local hot spot is located at to 0 < Z -5 
cm. which seems to indicate that film coolant is moved to other locations in 
the boundary layer by the vortex. 
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In Figure 36, the temperature field is presented for an expe- 
riment in which the heat transfer plate is heated to 40 2 C and film coolant to 
51 2 C. The temperature difference contours show the same qualitative 
trends as in Figure 35. In particular, the positions of many temperature 
gradients and the dimensions of the area affected by the vortex are about 
the same. Results in Figure 36 are different from those in 35 since : (1) it is 
not as easy to discern which part of the temperature field is due to film 
cooling only, (2) overall temperature differences are larger, and (3) local hot 
spots corresponding to film cooling locations are more easily discernable 
every 3 cm starting at Z = -12 cm. 

The temperature gradient at Z = -14 cm, corresponds to the 
spanwise edge of the film cooling for blowing rate of M - 0.98. 

c. Effect of spanwise vortex position. 

To study the effect of the vortex position relative to the film 
cooling injection, the vortex position was changed from A to B to C locations. 
Distributions of St/Sto are shown in Figures 33. 37 and 40, distributions of 
St/Stf in Figures 34, 38, 41 and distributions of T - Too are given in Figures 
35. 39 and 42. Regardless of the vortex position, the following phenomena 
were observed :(1) on the upwash side of the vortex a low heat transfer 
region exists, where St/Sto is less than for a film cooled turbulent boundary 
layer only, (2) on the downwash side of the vortex a wide region of high 
St/Sto exists, because of thinning of the boundary layer by the vortex, and 
(3) at X/d = 7.4, the film cooling dominates the flow and the vortex seems to 
have very little effect on the heat transfer. 
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Comparisons of Figures 33, 37 and 40, show that for low X/d 
values up to 17.5, a double peak in the St/Sto distributions is present for 
vortex positions A and B. At X/d = 33.6 and further downstream for vortex 
positions A and C, the double peak in St/Sto distributions starts to grow 
spanwise enlarging the portion of the plate covered by the hot spot. In these 
cases, St/Sto peaks are not as high as for the vortex position B case i Figure 
33 1 A maximum St/Sto of 1.05 was observed for the vortex' position b at 
X/d = 06.6, or about 30 percent greater than the St/Sto with film cooling 
only. Low St/Sto and St/Stf regions on the upwash side oi the vortex were 
observed for all three vortex positions. These were less significant for 
position C, while varying between 2- 6 percent for the other two vortex 
positions. 

St/Stf distributions from Figures 34, 38 and 40 show that the 
location, shape and size of the hot spot vary depending on the vortex span- 
wise position. In Figure 34, for vortex position B, this hot spot is located 
from the middle length of the plate, growing in size as the vortex develops 
downstream. For vortex position A, the hot spot is shifted upstream and is 
not as wide, nor as long as for position B. For vortex position C, the hot spot 
begins early upstream but persists further downstream, but it does not 
reach the end of the test section, and it is very narrow. A region of high heal 
transfer is observed in all three cases, but this area is larger for vortex posi- 
tion A, covering about 2/3 of the test plate. A heat transfer gradient is 
aligned along the vortex center for all three vortex positions, but it is 
steeper and more persistent for vortex position B. The region of low heat 
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transfer, where Si/Stf is less than 1, is largest for position B and smallest for 
position C. 

Figures 35, 39 and 42 show the T - Too temperature field in 
degrees C for X/d - 41.9. For these results the heat transfer plate was 
maintained at the freestream temperature and the coolant was heated to >1 
Q C, (unheated plate). Overall qualitative trends for all three vortex positions 
are similar: coolant is lifted away the wall by the upwash side of the vortex , 
on the downwash side of the vortex, secondary flows cause the freestream 
fluid to be located very near the wall, making the boundary layer very thin 
and minimizing the effect of the film cooling. A significant spanwise gradient 
in temperature exists which is located near the vortex core, and extends 
about 5 cm from the wall. As expected, the spanwise location of this gradi- 
ent changes with vortex position. The effect of the downwash side or the 
vortex is most significant for vortex position B. For this case Figure 35a 
shows that the boundary layer is thinned such that it only extends 2 cm 
from the wall. Significant changes also occur near the upwash side of the 
vortex where secondary flows from the vortex convect film coolant away 
from the wall. The most significant changes in this part of the flow occui for 
vortex position C. 

3. Free stream 1*> m/s. blowing ratio. M = 0.47. Effect of 

low blowing ratio. 

Measurements were made at Uco =15 m/s, with a blowing ratio 
M of 0.47, in order to determine the effect of low blowing ratio on heat 
transfer. These tests were also conducted to obtain additional verification of 
general conclusions from results obtained at 10 m/s and M = 0.98 blowing 
ratio. For these tests, vortex generator # 2, was used at positions A, B and C. 
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Resuits are shown in Figures 43 - 48. Distributions of St/Sto are shown m 
Figures 43. 45, and 47, and distributions of St/Stf are given in Figures 4 4 4b 
and 48. 

From a general qualitative point of view, the results at 15 m/s and 
M = 0.47 are similar to those of 10 m/s and M = 0.98. The most significant 
quantitative difference is that St/Sto maxima are higher, and St/Sto minima 
are lower for M = 0.47 In addition, the peaks of local heat transfer seem to 
increase more rapidly with dow r nstream distance. At the downstream end oi 
the plate where X/d = 96.6. local St/Sto for M ■•= 0.47 may be as high as 112 
compared to about 1.05 for M = 0.98. The rapid downstream growth of these 
maxima becomes more apparent when one considers the film cooled bound- 
ary layer alone, since St/Sto ratios at M = 0.47 are lower than those at M = 
0.98. 

At M = 0.47, the effect of changes in the spanwise position of the 
vortex are much less significant than at M = 0.98. At higher blowing ratio, 
many quantitative changes in the local heat transfer result as the vortex 
position is changed. In particular, primary and secondary peaks of St/Sto 
were observed at X/d = 17.5 and 33.6 for vortex position B, and X/d = 54 6 
and 75.6 for vortex position C. At the M = 0.47 blowing ratio, the shapes of 
St/Sto distributions for different vortex positions are very similar, especially 
for X/d = 54.6, 75.6 and 96.6. Generally, quantitative changes for M = 0 47 
occur only in the spanwise locations of the St/Sto maxima and minima. At 
X/d = 7.4, the effect of the vortex position changes is minimal and the 
spanwise heat transfer rates are mostly affected by the film cooling. At X/d 
= 17.5 both primary and secondary St/Sto peaks are present for all three 
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vortex positions. The only significant St/Sto changes with vortex position oc- 
cur at X/d - 33 6 where a small secondary St/Sto peak is observed when the 
vortex is located at position A only. Such behavior indicates that the effects 
of film cooling at M = 0.47 do not persist as far downstream as at higher 
blowing ratios. In addition, the vortex seems to more completely dominate 
the flow field as the blowing ratio decreases. 

4 . Free stream 10 m/s. blowing ratio. M 1,26. FiTect o f 

high blowing ratio. 

In order to study the effect of high blowing ratio on the film cooled 
turbulent boundary layer with embedded vortex, measurements were made 
at Uco = 10 m/s. with M = 1.26 and the vortex at position B. Results are 
shown in Figures 49- 52. 

The St/Sto and St/Stf data in 49 and 50 show the same overall 
qualitative trends shown by measurements at M = 0.98, 0.47, and 0.68 : 
augmented heat transfer near the downwash side of the vortex and reduced 
heat transfer near the upwash side. Also, at X/d = 7.4, film cooling rather 
than the vortex dominates the spanwise variation of heat transfer. The most 
important quantitative differences between these results and those at other 
blowing ratios are the magnitudes of St/Sto with film cooling only. Figure 26 
shows M = 1.26 St/Sto data at 10 m/s which is, on the whole, greater than 
data at M = 0.68, but just lower than results for M = 0.98. Results presented 
in Figure 49 reflect the same overall trends. 

Details of the spatial variations of normalized Stanton numbers in 
Figures 49 and 50 are now discussed. At X/d = 17.5, 33.6 and 54.6 double 
peaks in St/Sto distributions are present. These seem to result from inter- 
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action between cooling jets and the vortex, occurring atZ = -2 — 3 cm and at 
Z = 2 - 3 cm. 

Contours of T - Too (QC), at X/d = 41.9 in Figure 51 show that the 
peak near +2-3 era results from the downwash of free stream fluid to near 
wall regions by the vortex. The corresponding St/Sto peak increases m 
magnitude as the the boundary layer convects downstream reaching mag- 
nitudes as high as 1.05 at X/d 75.6 and 96.6. At these downstream locations, 
augmented St/Sto cover a large spanwise portion of the film cooled test sur- 
face extending from -5 cm to +5 cm (Figures 49 and 50 1. Magnitudes ot 
St/Sto maxima in these hot spots at M = 1.26 are about the same as those at 
M = 0.98. 

Figure 50 shows the heat transfer rate over the test plate in terms 
of St/Stf . Here, the hot spot is located from X/d = 54.6 to X/d = 96.6. This 
region is very large, confined to the downwash side of the vortex. The tem- 
perature gradient on the upwash side of the vortex is not very steep in some 
streamwise locations. A comparison between Figure 50 and Figure 32 where 
M = 0.68, show that, the hot spot is present at the same spanwise and 
streamwise location in both cases, but in 50 the hot spot is thinner, and the 
area of low St/Sto on the upwash side of the vortex is larger. 

Temperature contours at X/d = 41.9 in Figure 51, show that the 
St/Sto peak atZ = -2 — 3 cm. lies beneath a region of large temperature gra- 
dients very near the vortex center. Figure 52d. shows strong secondary 
flows at this location which apparently sweep coolant further from the vor- 
tex center (to smaller Z), such that it collects near -5 Z -8 cm Contours of 
vorticitv magnitude, calculated from these secondary flow vectors for M = 
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1.26 vortex position B, are shown in Figure 52e. These show that (1) the 
core of the primary vortex is located at Z = -2 cm., Y = 3 5 cm., and (2) that 
the secondary peak in heat transfer is not a result of any form of secondary 
vortex. Augmented mixing is initially created by the secondary flows 
beneath the vortex center, resulting in locally higher heat transfer. This se- 
condary peak shown in Figure 49 at X/d - 17.5, 33.6 and 54.6 eventually 
becomes indistinguishable at larger X/d equal to 75.6 and 96.6. The se- 
condary peak is bigger and more persistent than at M = 0.98, indicating that 
the shear layer interaction causing higher heat transfer is dependent upon 
the mass flux of film cooling. 

Another important feature of the results in Figure 5 1 is the gradient 
of heat transfer present at Z = -10- -12 cm. for 0 < Y < 5.0 cm. This gradient 
is present at the edge of the film cooled region : coolant was injected from 
only 7 holes extending across a span from -9 cm. to +9 cm. 

Additional flow field information at X = 1.48 m, X/d = 41.9 is given 
in Figures 52a.-c for 10 m/s, M = 1.26 and vortex position B. Data was taken 
using a five-hole pressure probe, [Ref. 1.: pp. 14-151, to measure pressure at 
800 points in a spanwise plane in order to further investigate the effect of 
the high blowing ratio in the turbulent boundary layer with embedded vor- 
tex. These results are consistent with those of Evans, lRef.1]. They show : ( 1 ) 
that high velocity, high total pressure fluid is swept near the wall on the 
downwash side of the vortex, (2) a velocity and total pressure deficits exist 
at the center of the vortex, and (3) low velocity fluid is swept away from the 
wall near the upwash side of the vortex. Vorticity contours in 52e. show a 
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large region of negative vorticity in the upwash region, (-8 < Z < -4 cm). In 
addition, two other regions of negative vorticity result from the secondary 
flows present. 
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IV. SUMMARY AND CONCLUSIONS. 



Heat transfer measurements were made of turbulent boundary layers 
with film cooling and embedded vortices for freestream velocities of 10 and 
15 m/s . To obtain the data, a heat transfer test surface was designed and 
developed to provide constant heat flux over its area, with 126 embedded 
thermocouples for detailed spanwise resolution of temperature. Extensive 
qualification tests show mat the surface gives excellent spatially resolved 
heat transfer coefficients over an area of 43.8 1 5 cm x 1 1 1 1.76 cm. 

Baseline measurements show excellent agreement with Stanton number 
correlations for a flat plate with constant wall heat flux and unheated start- 
ing length. Results of turbulent boundary layer with embedded vortex show 
excellent agreement with Joseph, [Ref. 3.1, and with the literature. Results of 
turbulent boundary layer with film cooling at different blowing ratios show 
expected trends, consistent with the data of Goldstein and Yoshida, (Ref. 1 8). 

Longitudinal vortices cause significant changes in heat transfer and 
structural characteristics of film cooled turbulent boundary layers. Heat 
transfer augmentations as large as 30 percent, and reductions as high as 10 
percent were observed to persist as many as 23 boundary layer thicknesses 
or 96 film cooling hole diameters downstream of film cooling injection loca- 
tions. The effects of the embedded vortex on heat transfer in film cooled 
boundary layers are significant and important : 

1.) Near the downwash side of the vortex the heat transfer is augmented, 
vortex effects totally dominate the flow behavior, and the effects of 
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the film cooling are negated. " Hot spots ” will exist for blowing ratios 
ranging from 0.47 to 1.26, and vortex circulation to freestream velocity 
ratios of about 1.6 cm. 

2. ) Near the upwash side of the vortex, the coolant is lifted off the wall and 

pushed to the side of the vortex, increasing the surface area protected 
by the film cooling. 

3. ) Changing the position of the vortex with respect to the film cooling jets 

results in significant local quantitative changes in heat transfer occur, 
even though the overall qualitative trends remain unchanged. 

4. ) Near film cooling injections locations for X/d up to 7.4, the film cooling 

dominates the flow behavior and the vortex seems to have very little 
effect on spanwise variations of heat transfer. 

Results (1) and (2) are mostly a consequence of the intense secondary 
flows produced in the plane perpendicular to the axis of mean vorticity. 
These results are consistent with those obtained by Joseph, [Ref. 3.: pp. 54], 
and Evans, [Ref. 1.: pp. 28]. 

At high blowing ratios: 

1. ) A double peak in the St/Sto distributions was observed to occur 

between X/d = 7.4 and 54.6. 

2. ) The change in spanwise position of the vortex in relation with the film 

cooling jets affects the magnitude, shape and spanwise position of 
St/Sto peaks, and 

3. ) Secondary St/Sto peaks become higher in magnitude and more 

persistent with respect to downstream development as the blowing 
ratio increases. The double peaks observed at M = 1.26 were not due 
to a secondary vortex, but to an interaction between the vortex and 
film cooling which depends on blowing ratio. 
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At low blowing ratios: 

1. ) The change in spanwise position of the vortex has very little effect on 

local heat transfer distributions except to change the locations of 

Stanton number minima and maxima, and 

2. ) St/Sto distributions exhibit only one peak which increases in magni- 

tude with downstream development. 

It is recommended that flow visualization study of the interaction of the 
vortex and film cooling be conducted in order to enhance the understanding 
of some of the complex phenomena observed during the course of this 
study. 
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Temperature Bounadary Layer 
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VORTEX GENERATOR *2 POSITION 
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Fig. 2. Vortex Generator Position. 
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Fig. 3b. Detail of the Injection Plenum and Vortex Generator 
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CD VS. REYNOLDS NUMBER 
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REYNOLDS NUMBER *E + 3 

Discharge Coefficient vs. Reynolds Number 




Fig. 5a. Heat Transfer Test Plate 
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Fig. 5b. Detail Of Heat Transfer Plate 
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All dimensions in cm. 

Fig. 6. Thermocouple Placement on Plate 
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THERMOCOUPLES (.762 mm) 




Fig. 7. Cross Section of Test Surface 
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Temp Plate - Temp. Amt) ( ? C) 

Fig. 8. Conduction Losses 
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Z (CM) 

FREE STRERM 10 M/S, NO FILM COOLING 

Spanwise Variation of Heat Transfer Coefficients. 
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REx * 10^6 

FREESTRERM 10 M/S, NO FILM COOLING 

Spanwise Averaged Stanton Numbers. 
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11. BOUNDARY LAYER, 10 M/S NO FILM COOLING 
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Z (CM) 

FREE STRERM 15 M/S, NO FILM COOLING 

Spanwise Variation of Heat Transfer Coefficients. 



AVERAGED STANTON NUMBERS 

DRTE - 7 1?87 . 1702 



Q 

« 

r> 




csneoneoneotneoneanco 
ior>r>'c^cnoc\jro~*-* Q 



£-v01* ON NOlNblS 






75 



REx * 10^6 

13 . FREESTRERM 15 M/S, NO FILM COOLING 
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Fig. 14. FREE STRERM !0 M/S, NO FILM COOLING 
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rig. 15. FREE STREAM 15 M/S , NO FILM COOLING 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 
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Fig. 16. FREE STRERM 10 M/S, WITH FILM COOLING 
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STRNTON NUMBER RRTIOS 

DATE - 725B7 . 1 4 
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Fig. 18. FREE STRERM 10 M/S, WITH FILM COOLING 
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Fig 19 FREE STRERM 15 M/S, WITH FILM COOLING 
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Fig. 20. FREE STRERM 15 M/S, WITH FILM COOLING 
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REYNOLDS No* 10' 

FREE STREAM 10 M/S 
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REYNOLDS No* 10' 

FREE STREAM 10 M/S 



FILM COOLING 
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REYNOLDS No* 10' 

FREE STREAM IS M/S 
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Fig. 27. Free Stream Velocity IS m/s 
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Fig. 28. Free Stream Velocity 10 m/s 
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TEMPERATURE PROEILES 
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STRNTON NUMBER RATIOS 

DATE - 7280? . 1742 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79cm 
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FREE STRERM 10 M/S. WITH FILM COOLING, VORTEX GEN. RT 4.79cm 
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STANTON NUMBER RATIOS 

DATE - 72007 . 1742 
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STANTON NUMBER RATIOS 
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FREE STRERM 10 M/S, WITH FILM COOLING, VORTEX GEN. RT 4.79cm 

Fig. 31e. Spanwise Variation of Stanton Number Ratios. 
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SURFACE CONTOURS 

DOTE - 72807 . 1742 
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FREE STRERM 10 M.'S, WITH FILM COOLING, VORTEX GEN. #2 FIT 4.79 cm 

Fig. 32. Surface Contours. 
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STHNTON NUMBER RATIOS 

DATE - 80987.1943 




101 



-20.00 -19.00 -10.00 -S.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STRERM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 RT 4.79 cm 

Fig. 33b. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 13.00 - 10.00 - S .00 0.00 3.00 10.00 13.00 20.00 
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FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 4.79 cm 

Fig. 3oC. Spanwise Variations of Stanton Number Ratios. 
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-20. 00 -IS. 00 -10.00 -5.00 0.00 3.00 10.00 IS. 00 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 4.79 cm 

Fig. 33d. Spanwise Variations of Stanton Number Ratios. 
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FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 4.79 cm 

Fig. 33e. Spanwise Variations of Stanton Number Ratios. 
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-20.00 -IS. 00 -10.00 -S.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 4.79 cm 

Fig. 33f. Spanwise Variations of Stanton Number Ratios. 
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FREE STRERM 10 M/S, WITH FILM COOLING, VORTEX GEN. #2 RT 
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TEMPERATURE PROFILES 
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Fig. 35d . EMBEDDED VORTEX t0 M/S WITH FILM COOLING 
X«2.4B0 M, UNHEflTED PLATE .VORTEX GEN. AT 
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STRNTON NUMBER RATIOS 

DATE - 00487 . 1 529 
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Z (CM) 

FREE STRERIi 10 M/S, NITH FILM COOLING, VORTEX GEN. *2 AT 3.52 cm 

Fig. 37a. Spanwise Variations of Stanton Number Ratios. 
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— 20.00 - 19.00 - 10.00 - 9.00 0.00 9.00 10.00 19.00 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. #2 AT 3.52 cm 

Fig. 37b. Spanwise Variations of Stanton Number Ratios. 
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FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 3.52 cm 

Fig. 37c. Spanwise Variations of Stanton Number Ratios. 
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FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 3.52 cm 

Fig. 37d . Spanwise Variations of Stanton Number Ratios. 
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STANTON NUMBER RATIOS 

DATE - 80487 .1529 
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** 20.00 - 19.00 - 10.00 - 9.00 0.00 9.00 10.00 19.00 20.00 

Z (CM) 

FREE STRERM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 RT 3.52 cm 

Fig.37f. Spanwise Variations of Stanton Number Ratios. 
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FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. #2 AT 3.52 cm 

Fig. 38. Surface Contour. 
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STRNTON NUMBER RRTIOS 

DRTE - 00407.1553 
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- 20.ee - 19.00 - 10.00 - 9.00 0.00 9.00 10.00 19.00 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT G.06 cm 

Fig. 40a. Spanwise Variations of Stanton Number Ratios. 



STRNTON NUMBER RRTIOS 

DATE - 00487 . 1 553 
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STANTON NUMBER RATIOS 

DATE - 80487 . 1 553 
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-2e.ee -i3.ee -ie.ee -s.ee e.ee s.ee ie.ee is.ee 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 6.06 cm 

Fig. 40c. Spanwise Variations of Stanton Number Ratios. 
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- 20 .ee -t 5 .ee -ie.ee - 5 .ee e.ee 5.00 ie.ee is.ee 20.00 

2 (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 6.06 cm 

Fig. 40d. Spanwise Variations of Stanton Number Ratios. 



STRNTON NUMBER RRTIOS 
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- 20.00 - 19.00 - 10.00 - 9.00 0.00 9.00 10.00 19.00 20.00 

Z (CM) 

FREE STRERM 10 M/S, WITH FILM COOLING, VORTEX GEN . *2 RT 6.06 cm 

Fig. 40e. Spanwise Variations of Stanton Number Ratios. 
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-2e.ee -is.ee -ie.ee -s.ee e.ee s.ee ie.ee is.ee 2e.ee 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. *2 AT 6.06 cm 

Fig. 40f. Spanwise Variations of Stanton Number Ratios. 
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Surface Contour. 
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Z (CM) 

42 . EMBEDDED VORTEX 10 M/S WITH FILM COOLING 
X- 1.4 00 M, UNHERTED PLATE, VOTEX GEN. AT 6.06 cm 




STRNTON NUMBER RRTIOS 

DATE - 72907 . 1308 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79cm 

Fig. 43a. Spanwise Variations of Stanton Number Ratios. 
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STRNTON NUMBER RRTIOS 

DRTE - 7298? . 1308 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79cm 

Fig. 43c. Spanwise Variations of Stanton Number Ratios. 



STANTON NUMBER RATIOS 

DATE - 72987. 130B 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79cm 

Fig. 43d. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79cm 

Fig. 43e. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79 cm 

Fig,. 43f. Spanwise Variations of Stanton Number Ratios. 
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FREE STRERM 15 M/S, WITH FILM COOLING, VORTEX GEN. #2 RT 

Fig. 44. Surface Contour. 



STRNTON NUMBER RATIOS 

DOTE - 7290? . 1345 




137 



- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STRERM 15 M/S, WITH FILM COOLING, VORTEX GEN. RT 3.52 cm 

Fig. 45a. Spanwise Variations of Stanton Number Ratios. 



STANTON NUMBER RATIOS 

DATE - 72967 . 1345 
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-2e.ee -is.ee -ie.ee - 5 .ee e.ee 5.00 ie.ee 15.ee 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 3.52 cm 

Fig. 45b. Spanwise Variations of Stanton Number Ratios. 
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-20 . 00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 3.52 cm 

Fig. 45c. Sparwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 3.52 cm 

Fig. 4 Sd . Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 3.52 cm 

Fig. 45e. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 3.52 cm 

Fig. 45f. Spanwise Variations of Stanton Number Ratios. 
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FREE STRERM 15 M/S, WITH FILM COOLING, VORTEX GEN. #2 RT 

Fig. 46. Surface Contour. 



STANTON NUMBER RATIOS 

DflTE - 7298? .1415 




S3 

S3 

• 

(9 

OJ 



(9 

<9 

in 



Q 

(9 

• 

S3 

(9 

(9 

in 

(9 

S3 

• 

o 



C9 

S3 

• 

in 

i 



(9 

S3 

• 

S3 

T 



S3 

S3 

in 

I 



S3 

S3 

d 

OJ 

I 



6 

u 

U3 

CS3 

CD 



Z 

u 

(J 

X 

u 

H- 

q: 

o 

> 



‘Si 

o 



03 

U 



?-■ 

o 



6 

3 



z 

u 



N 



(3 

z 

l-H 

_i 

o 

0 

u 

z 

_] 

1-4 

L. 

1 



G 

O 

4-> 

c 

03 

<p 

CO 



c 

o 



03 

u 

03 

> 



^ O 

cn 2 

i * 

2- C 

03 

in g* 

-4 CO 



Z 

CE 

U 

QL 

K 

cn 

u 

u 

a: 

u. 



03 

nr 

&0 

•H 

Uh 



144 



STRNTON NUMBER RATIOS 

DATE - 72987 . M 15 
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STRNT ON NUMBER RRTIOS 

DRTE - 72987 .1415 
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- 20. 80 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 6.06 cm 

Fig. 47c. Spanwise Variations of Stanton Number Ratios. 



STANTON NUMBER RATIOS 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT B . 06 cm 

Fig. 47d. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 6.06 cm 

Fig. 47d. Spanwise Variations of Stanton Number Ratios. 



STANTON NUMBER RRTIOS 

DATE “ 7290? .1415 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 15 M/S, WITH FILM COOLING, VORTEX GEN. AT 6 . 0G cm 

Fig. 47f. Spanwise Variations of Stanton Number Ratios. 
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FREE STRERM 15 M^S , WITH FILM COOLING, VORTEX GEN. #2 RT 

Fie. 48. Surface Contours 



STANTON NUMBER RATIOS 

DATE - 00107. 1 54 1 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79 cm 

Fig. 49a. Spanwise Variations of Stanton Number Ratios. 



TRNTON NUMBER RRTIOS 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79 cm 



STRNTON NUMBER RRTIOS 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79 cm 

Fig. 49c. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.79 cm 

Pig. 49d. Spanwise Variations of Stanton Number Ratios. 
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- 20.00 - 15.00 - 10.00 - 5.00 0.00 5.00 10.00 15.00 20.00 

Z (CM) 

FREE STREAM 10 M/S, WITH FILM COOLING, VORTEX GEN. AT 4.73 cm 

Fig. 49e. Spanwise Variations of Stanton Number Ratios. 



STANTON NUMBER RATIO 
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-2Q.00 -13.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00 

2 (CM) 

FREE STRERM 10 M/S, WITH FILM COOLING, VORTEX GEN. RT 4.7 3 cm 

Fig. 49f. Spanwise Variations of Stanton Number Ratios. 
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Surface Contours. 



TEMPERRTURE PROFILES 

DATE- 8120?. 1730 
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STRERMWISE VELOCITY 
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EMBEDDED VORTEX 10 M/S WITH FILM COOLING , M“ 1.26, VORTEX GEN RT 

Fig. 52a. Streamwise Velocity Contour 
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l-'ig. 52 b. Total Velocity Contour. 
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EMBEDDED VORTEX 10 M/S WITH FILM COOLING, 
VORTEX GEN. RT 4.79 cm 



SECONDRRY FLOW VECTORS 
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EMBEDDED VORTEX 10 M/S WITH FILM COOL ING , M- 1 . 2G , VORTEX GEN. RT 4.79 cm 



STREflMW I SE VORTICITY 
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EMBEDDED VORTEX 10 M/S, VORTEX GEN. AT 
FILM COOLING, M- 1 .26 



APPENDIX B 



UNCERTAINTY ANALYSIS 




Kline and McGintock l Ref. I9|. This is the root sum square' method, where 
the uncertainty, 6 f . of some function F, is a function of the independent 
variables X a , according to the following expression: 



All uncertainties are based on 95% confidence Ieveis. To calcinate the 
uncertainty of the Stanton number, it is necessary to calculate first the 
uncertainty in the heat transfer coefficient, h. The following independent 
variable uncertainties were determined : 8 q = t 16 W/m 2 . 5j'v = t 0.5 Q C Spv 
= t 0. i °C. Here, the uncertainty of neat ioss by convection was basea on a 
5 % error in radiation losses and 2.5% error in conduction losses. The 
uncertainty of T w is higher than Too due to higher uncertainty in the 
calculation of the contact resistance. From these parameters the uncertainty 
of h was determined to be 4.8% or approximately ± 1.7 W/m 2 -(!. based on an 
h value of 35 W/m 2a C. 

To calculate the uncertainty of the Stanton number, besides the uncer- 
tainty of the heat transfer coefficient, the following independent variable 
uncertainties w r ere determine : 8pco = i 0.0 i Kq/rrP, 5(jo© = r0.5 m/s and Cp = 



( n 



3 




leqn. CP 
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r ! 0 J/Kc C 'K. The uncertainty of C p was based on the assumption of constant 
properties. From these parameters the uncertainty of the Stanton namner 
was calculated to be t 6.04%, or approximately 1.87 x l O' 4 based on a typical 
Stanton number of 0.0035 The uncertainty of St/Sto is derived from the 
parameters already mentioned and was estimated to be ± “77 or t 0 0“ 
based on a St/Sto of 0.75. 



To calculate the uncertainty of the blowing ratio. VI. along the uncer- 
tainty values determined for Uoo and pco. the following independent 
variable uncertainties were: 5in n j - - 1.0 m/s, and §oi n j - - .02 'y rr . The 
uncertainty of U m j is higher than Uoo due to higher uncertainty of the 
rotometer and the area of the injection holes, and the uncertainty of pi n j is 
also higher than p=c due to higher uncertainty in the estimation of T m , from 
Tp. From these parameters the uncertainly of M was calculated to he - 8 d-t 
or about i 0.087c for a M value of 0.98. 
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APPENDIX C 



SOFTWARE 



The following programs are listed 



-t't-, r> 

'« CA\0 



cfltffiV OaiUflCd £$UlYlduOu i Of CULluUCUofl iUSScS 



ST 



\ vtr r i 

n..\fL 1 



heal transfer prc"**' 
readings from DAS., 
files in floppy disk 



m 10 acquire thermocouple 
and store information in data 



STANFC 2 : heat transfer program to calculate Stanton numbers 
for ail flow conditions 



ACQTPRO : heat transfer program to acquire temperatures from 
automatic traversing device for temperature surveys 

PLOTRUN : program to read a data file to plot temperature 
contours 



PTSTAV : programs read a data file to plot spanwise averaged 
Stanton numbers vs. Reynolds number 

PTSTLC : reads a data file to plot spanwise local heat transfer 
coefficients 

PLSTRTIO : plots spanwise St/Sto for film cooling only. 

PLSTRVOR: plots spanwise variations of St/Sto for embedded 
vortex data only 

PLSTRVY. plots spanwise variations of St/Sto for film cooling 

only and St/Sto for film cooling and embedded vortex 
by rows. 



166 



Sl’RFCONT: plots surface contours of St/Stf 



PLSTRFC 



plots s panwise variations of local St/Sto for film 
cooling data only 



Name of variables used are intended to be self explanatory 
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IQ 
20 
30 
40 
53 
SG 
70 
SG 
30 
133 
1 13 
I Z9 
i j3 
; 40 
150 
1 5G 
IS I 
\ s z 

1 73 
180 
130 
200 
Z I 0 
ZZG 
230 
Z 40 
Z5G 
263 
Z73 
283 

2 3G 
300 
310 

T r\ 
<J £~ ’Ll 

330 
343 
350 
38G 
373 
333 
333 
403 
413 
423 
430 
443 
453 
4a 'Ll 
470 
4SB 
480 
530 



i PROGRAM ENERGc 
I 

! THIS PROGRAM ACQUIRES MULTIPLE CHANEL THERMOCOUPLE DATA 
! AND PERFORMS ENERGY BALANCE TO ESTIMATE CONDUCTION LOSSES 
•LI5RANI/0RTIZ VERSION, JUNE ! SSS 
! 

DIM E( 233 ) ,T( 200 ) 

» 

! CHANNELS 0-73, i GO- i 46 
! COFFER CONSTANT AN THERMOCOUPLES 



PRINT 

PRINT "ENERGY BALANCE RESULTS" 

PPI,MT " 

! ENTER POWER IN (WATTS) 

DISP “( HIT .(CONTINUE'-" >“ 

~ rt \J Z ul 

2 ISP "ENTER SUPPLY CUPP ENT IN AMPS, AND VOLTA.SE IN VOLTS' 
INPUT nr^p 5 ,Vo I is 
PR I NT 

PRINT "CURRENT = *' ?Ar*ps , “ VOLTS - ’Waits 
PRINT 

PRINT "TEMPERATURE RESULTS'* 

P R I N T 

^ave— 3 . 

PRINT 

PGR 1 = 1 TG 1 3 
OUTPUT 70S ; “A I " ; I ; “ VT 1 11 
ENTER 709 ; X 
EC I > =X * 1 333300 . 

TC I /«.01SZ05+.02S848*E( I >- .000000531 *El I )*E( I > 
Tave*Tsve+Tv I ) 

PRINT USING 33G;I,£CI JJCI) 

IMAGE DDD ,3X f DDDDG.O ,3X ,0000,00 
NEXT I 

PRINT " + 

PRINT 

PRINT 

FOR 1=43 TO 5S 

OUTPUT 70S 5 " nl “ } I ; " VT1 " 

ENTER 70S : X 
E< I ) = X * i 00 0 000 . 

T< I )= . Q rS20Sr . 32E34S *Ev I >-.00000053! *E< I 5*E< I 5 
Tave=Tave+T( I > 

PRINT USING 453 ; I ,E( I > ,T( I > 

IMAGE DDD.3X .DDDDD.D ,3X .DDDD.DD 
NEXT I 

PRINT " **t + *<****^/*»*^*** + »4^*+» ) f^**t*r + ***** + *t++*+t* N 

PRINT 

Tave=Tave/4S . 

PRINT "AVERAGE PLATE TEMP. MEASURED (DEG C> »";Tave 
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540 PRINT 

550 ! GET T AMBIENT 

SS0 OUTPUT 70S ; "AI “ ; 0 ; " VT ! “ 

570 ENTER 709 ;X 
580 E<0>«X*1000 

530 T<3>=26. 573*5(0)-! . 936879*E< 3 >*E( 3 ) + . 9S795*E< 3 >=( 0 )*E< O )- . 25 1 27?*E( 0 >*E( 0 

)*£( 3 )•♦£( 8 ) 



613 


Tarcb=T C 0 ) 




620 


PRINT "AMBIENT TEMP 


. = " ; Tawb 


630 


PRINT 




S40 


Tdi f f=Tave-Tanb 




553 


PRINT "TPLATE-TAKB 


= 4 ; Tdi f f 


663 


PRINT 




670 


! GET AVERAGE TEMP I 


N INSULATION 


650 


Tina 1 =3 . 




S93 


FOR 1=147 TO US 




703 


OUTPUT 733 : " AI " ; I ; " 


VT 1 " 


713 


ENTER 70S i X 




Tin 

/ i. 


E< I )=X*1003. 




730 


T( I >=26 . 573 *E( I 5-1 . 


S36379*EC I > r E 


> *E\ 


I > * E ( I > 




750 


Ti ns 1 =T i ns 1 +T< I ) 




763 


! PRINT USING 463; I, 


E( I > ,7 < I ) 


770 


NEXT I 




733 


Tins ' =T ins 1 /3 . 





■ e; i 



1 



790 T m52=0 . 

503 FOR 1=153 TO 152 

310 GUTFUT 703; "AI " ; I ; “UT! " 

S23 ENTER 705 ;X 
833 E<IJ-X*1000. 

340 T ( I >=26 . 573*E< I >-l .33G8?9*Ev I >*E< I )-,9378S*E( I > 
>*E( I )*E( I ) 

850 (PRINT USING 450 ; I ,E< I > ,T< I) 

370 Tins2=Tins2+T\ I ) 

883 NEXT I 

890 Tin52 = Ti ns2/'3 . 

300 PRINT 

S 1 3 PRINT "TINS 1 = " ; Tins 1 , " TINSZ «’';Tins2 
323 PRINT 



E< I )«E< I >-.261 Z77+E; I >*E< I 



930 1 ENTER VALUES OF THERMALCONOUCT . , AREA , AND THICKNESS OF INSULATION 

940 K= .04 ! W/M DEG C 

353 Q:<= . 0254 !M 

960 A= . 4337 ! M " 2 

973 ! CALCULATE HEAT FLUX THROUGH INSULATION 

3S3 Qins=K *A*< Tins i -Tina2 )/Dx 

330 ! SET 0 IN ( WATTS 5 

1000 Qin=VoI i 5 ♦Amps 

1010 Qcond=Qi n-Qins 

1020 PRINT “Q IN = " ; G i r. , " QINSULATION »’?Qins," QCONDUCTIGN =";Qcond 
1033 PRINT “ * * ■* *************** * *** ^ ***♦■***»■*■***•*•-*♦* 4 ». * + »• * '■ 

1040 STOP 



1053 END 
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13 


i pros; 


23 


i 


33 


t THIS i 


40 


fCREAT! 


S3 


1UFDATI 


S3 


» 


73 


( 


3 3 


MASS 3 


90 


CREATE 


133 


A 3 S 1 9N 


no 


CREATE 


1 23 


ASS IGN 


1 30 


DIM E ' 


140 


FRIKTEi 

l 


! S 0 


PRINT 


1 73 


INPUT 


130 


IF An 5 


i r*> .-i 

1 C ‘J 
7 ftp, 


• 


213 


• CHANE 1 

icoppe; 


223 


j 


~7 m {* 
u» J «J 


PRINT 


2*0 


INPUT 


25*0 


PR [NT 


2 S3 


PR I NT 

j 


2 SO 


i ENTER 


250 


PRINT 


300 


INPUT 


313 


PRINT 


320 


INPUT ! 




PRINT 


343 


PRINT 


352 


i ENTER 


3 SO 


PRINT 


370 


INPUT. 


JGU 


PRINT 


333 


INPUT 


430 


PRINT 


410 


INPUT 


420 


PRINT 


433 


PRINT 


440 


! 


4S3 


OI3P “ 


4S0 


PAUSE 


473 


PRINTE 


480 


PRINT 


450 


PRINT 


533 


PRINT 


510 


PRINT 


523 


PRINT 



MASS STORAGE IS " ; INTERNAL ,4 , ' 
CREATE BOAT “ T0AFC2S " ,252 ,S 
ISN OPa t i TO “ TGAFC2S 
BOAT * I DAF C2S '' , I I ,3 

gp^i h_: to , ‘ida c c:e“ 

30 ; .Tv 2 S 0 5 



I OF EM FILE TO RETRIEVE BASELINE TEMrERATU 



OPEN FILE TO RETRIEVE INPUT CAT; 



! SUBROUTINE FOR FILMCCOLING DATA 



v '> • r- r 

;3~ i wt 



“RUN S < MONTH, DAY, YEAR. HOUR .MINUTES) -"jRunnc 



)=" ;Ta«b 



-♦-f-***-***-* 
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530 PRINT "Mo T DEG G" 

54-0 ! ACQUIRE THERMOCOUPLE READINGS 

553 PRINT 

550 FOR 1=] TO 79 

570 OUTPUT 70S;"AI"iIi*UT1" 

5S0 ENTER 709 ?X 

530 El I >=X * i 300300 

S00 T < I >-.Q1S20S+.0Z5S4S*E{ I >-.000000531 *E< I 3 -E\ I > 

E i 0 PRINT USING 520? I ,T( I ) 

S2G IMAGE ODD ,4X ,DQ.CQ 



533 


OUTPUT ®Fath_i 


540 


NEXT I 


553 


! 


553 


FOR 1=] 30 TO I 


/*• •= * r\ 
□ i 


OUTPUT ?2S:”AI 


SS0 


ENTER 789 iX 


S S3 


E\ I ) = X- : 303203 


700 


T( I /= . G I 325+ . 3! 


710 


11=1-23 


720 


PRINT USING S2i 


750 


OUTPUT @Path_l 


740 


NEXT I 


753 


! SET AMBIENT Tl 


76 S 


OUTPUT 733; "AI 


t ( V 


ENTER 733 ; X 


7SG 


E( 3 >=X » 1 333 . 


73*3 


T ( 8 )=2S . 573^Ev 1 


633 


Tcharib = T v 3 ) 


3 ! 0 


PRINT “ THERMQC 1 




PRINT 


830 


! GET THE REST i 


S40 


FOR 1 = 147 TO I 1 


350 


OUTPUT 70S; "A I 


850 


ENTER 733 ; X 


870 


E< I )=X * ! 303 , 


883 


T ( I )=2S . 573*-E( 


330 


PRINT USING 52' 


930 


OUTPUT QPath_1 


31 3 


NEXT I 


920 


ASSIGN OFat h_i 



s " ; i charcb 



330 
343 
350 
35 S 

•-* r\ 

S(W 

930- 

930 

1000 

J010 

I 

1030 
1 v543 

jasa 
1 680 



! TRANSFORM TO SI UNITS 
| 

Parcb-Farch*33£5 . 82 
Dei tap-Qel tap*Z48 .7 
F 5 tanp^Ix J 47 >+2?3. IS 
Ro“Panb/( 237-Fs i enp ) 

Uinv={ 2 ►Dei tap/Ro )' ,5 
|. . 

PRINT 

OUTPUT OF a t h_2 ; Runnc f ftmp $ /Jolts ,Parnb , Del tap , Tamb .Tchamb f Ui nf 3 Rg , Fs t erop , An 
AS S I 3h @Pa t h_2 TO 4 

f - 

GOTO I 1 90 

f SUBROUTINE TO CALCULATE FILM COOLING DATA 



i pascals*: N/r-rz ) 

! PASCALS* N/M~2 ) 

! DEG KELVIN 
f AIR DENSITY (KG/M3) 

1 FREE STREAM VELOCITY <M/S> 
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1070 PRINTER 13 1 

1030 PRINT "ENTER FLOW % FROM ROTOMETER " 

1G30 INPUT Fiw 

11.00 PRINT "ENTER PLENUM FRESSURE -DIFFERENCE (IN H20)“ 

11,10 INPUT Oelpi 

1122 Delpi=Deipi*243.7 

1130 Flw=Flw* .000S3345E 

1140 CREATE BOAT "FILQT26 ’,2,3 

1153 ASSIGN @Patn3 TO "PILOTS^ " ! OPEN FILE TO RETRIEVE FILNCOOLIN6 DATA 

11S0 OUTPUT SFaihPjFIw.OpIpi 
1170 ASSIGN @Path3 TO * 

1 130 GOTO 1 33 

1130 MASS STORAGE IS INTERNAL ,4“ 

1200 END 



( 
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10 
20 
30 
40 
50 
S0 
70 
80 
90 
100 
110 
120 
1 30 
U0 
150 

iee 

170 

188 
t r* r% 

* JkJ 

n (? r> 
L. !tX‘ 

210 

^ ^ 0 

233 

240 

25£ 



! FROGRAM STANFC2 

! 

IIQ M/S FREE STREAM VELOCITY , 13 INJECTION HOLES OFEN 
(THIS PROGRAM ACQUIRES MULTIPLE CHANEL THERMOCOUPLE DATA 
! CALCULATES HEAT TRANSFER COEFFICIENTS AND STANTON NUMBERS RATIOS 
! AND- FILMCOOLING PARAMETERS. 

! LIGRA.NI/ORTIZ/ JOSEPH VERSION , NOVEMBER 1S3G 
! UPDATED BY ORTIZ, JULY 1987. 



MASS STORAGE IS INTERNAL ,4 .. I" 

CREATE BOAT “STRFCVi 0" ,504 ,8 I OPEN FILE TO RETRIEVE STANTON Nos RATIOS 

ASSIGN SFath 4 TO “STRFCVI 0" 



DIM T ( 200 > 
DIM Ho (230 



,H( 200 ) , St ( 200 ) , X 1 ( S ) ,Hi (8 / , 3 1 1 (S' , X 2 ( 23 G ) , Z ( 203 5 . 3 1 f ( 
> ,Stc( 2O0 ) , X 1 a ( 3 ) ,X2a< 200 ) ,Za(2 00 ) ,Str(2C0) ,St fr(203 ) 



i 



PRINTER IS ! 



PRINT “ARE YOU USING VORTEX GENERATOR? (3.Y OR 4,N) 
INPUT Y 
IF Y-3 THEM 

PRINT “ENTER TYPE OF VORTEX GENERATOR " 

INPUT Uort 

PRINT “ENTER LOCATION FROM CENTER LIME IN c,V 
INPUT Loc 



2B0 
270 
280 
230 
300 
31 S 
323 
330 
340 
350 
350 
370 
330 
393 
430 
410 
420 
430 
443 
450 
480 
470 
480 
433 
500 
510 
520 



ELSE 

GOTO 290 
END IF 

PRINTER IS 701 

PRINT CHRS< 27 ) j CHR $( 33 ) ; CHR3( 1 33 > :CHR$< 53 > ; CHR2( 52 ) ; CHR'3( 70 ) 

PRINT CHRS( 27 ) jCHRSt 38 ) }CHR3< 108 ) ;CHR3( 49 >;CHRS( 75 1 

I 

ASSIGN SPathi TO “TDAFCIO" - .. • 

PGR 1=1 TO 133 
ENTER ©Path! *, T ( I > 

NEXT I 

ASSIGN aPathZ-TO “IDAFO10“ 

ENTER @Path2 sRunno ,Atnps ,Val ts ,Fanb ,Cal tap ,Tamb , Tcharcb .Uinf ,Ro ,Fstai"p ,Ans 
PRINT 

l 



DISP "(HIT < CONT I NUE > ) “ 

PAUSE 

IF Ans= i • THEN. 

PRINT "STANTON NUMBER RESULTS WITH FILM COOLING" 

else' 

PRINT "STANTON NUMBER RESULTS” 

END- IF • 

PRINT “**’**ir*#fr4***1***##')»1f***4#**********lHt***<f***<‘-********'*»* 
PRINT 

PRINT “RUN If (MONTH ,DAY .YEAR. HOUR , MIN >=" iRunno 
PRINT 

IF Y-3 THEN 
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S30 


PRINT "VORTEX GENERATOR IS No «Vort , “ LOCATED AT 


( on ) - !l ; Loc 


S43 


PRINT 




533 


ELSE 




530 


GOTO 3S0 




573 


END IF 




533 


IF Ans= t THEN GOTO 2243 




530 


! PRINT OUT DATA 




300 


PRINT "AMPS = u iAnp5,“ VOLTS -‘s Volts 




3 ! 3 


PRINT 




323 


PRINT "AMBIENT PRESSURE (N/M"2> = ";Panb," DELTA P 


(N/fT 2 ) i Dai tap 


533 


PRINT 




3 43 


! 




330 


PRINT “AIR DENSITY ( KG. -'M3 ) -*jRo,“ VELOCITY (M/3 1 


- " ; U i n f 


333 


PRINT 




573 


PRINT "FREE STREAM TEMP ,( DEG K " iFatspip 




SS3 


P R I N T 




533 


PRINT “THERMOMETER AMBIENT TEMPERATURE (DEG C) 


i afab 


733 


PRINT 




713 


! CALCULATE THE AVERAGE PLATE TEMPERATURE 




720 


| 




723 


T _ . _ 

i sj v e — tu . 




743 


FOR I •- 1 70. 125 




75G* 


T a v e ~ ‘ dveri [ I ) 




753 


NEXT I 




773 


T a v a = 1 a v c / i 2 S 




( C v 


PRINT "AVERAGE PLATE TEMPERATURE. MEASURED ( DEG C) 


iTdve 


i Z 


PRINT 




330 


!■ ENERGY cALANCh 




• /“» 
O l V 


| 




320 


T d i r r = f avc-To he fab 




S30 


Qcond=.oS3-r.3E4*Tdiff-.0jS*Tdiff /, 2 




S4G 


Tdvab3 3= Ta\c;T'27j. 1 5 




S50 


Tchaba-Tchanb+273 . 1 5 




<-» <-> 
oat? 


Qr ad= . QC0GGG02 1 E 5 ^ < Ta^/abs A 4-Tchabs ,N 4 ) 





-S70 Qi n=n(*ipa H'ol t a 
330 Ocsnv-Qin-Qcond-Qrad 
350 Dei t=Tavaba-Fstemp 
300 i 



PRINT “T PLATE - T P STREAM (DEC- C )=“ s 0= i t , "POWER IN ( WATTS )=" i Gin 
PRINT 

PRINT ‘CONDUCTION LC33 (WATTS)-* iQcend ."CONVECTION LOSS (WATTS >-* {Qconv 
PRINT 

PRINT "RADIATION LOSS (WATTS) =" ;Crad 
PRINT 

I WALL TEMPERATURES CORRECTIONS 



310 
323 
933 
540 
S53 
SS0 
573 
363 ! 

330 Cp=1 033 . 
? 330 Cr= .014 



i SPECIFIC NEAT FOR AIR 
! CONTACT RESISTANCE CORRECTION FACTOR 



1013 Dt»Cr*Qconv ! CELT A TEMP. DUE- TO CONTACT RESISTANCE 

1323 Q-Gconv/ . 4337 i CONVECTION LOSS CORRECTION 

1333 i 

1343 Tcavfcs=Tavabs-Ot 



1350 Del tc=Tca vbs-Fa t ewp _ 

1 BSC Theta=( T i jaba-Fstetip )/Del to 
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IS 10 

1 623 

1 530 

1640 

1 650 

1650 

1 670 

1630 

1630 

1700 

1710 

1720 

1730 

1 743 

1 7S3 

1763 

1770 

1 73G 

1750 

1300' 

1310 

1 3Z3 

1 833 

1340 



I860 

1870 



Z( M )-Z ( 1 ) 

X2(M) = X3 

FOR 1=11+1 TQ J *2 1 
X2( I )=X3 

Z( I >=Z( I-i >+1 .27 
NEXT I 
X3=X3+ . 2 
NEXT J 

FDR 1 = 1 TC 1 26 

OUTPUT @Fatn_4;Str< I ) ,5t fr< I ) ,X2< I J r Z C I ) 

PRINT USING !720d ,X2( I > , Z ( I ) , H < I ),3tr< I > ,5t fr( I ) 

IMAGE DOD ,2X ,00.00 ,3X ,3D,D0 ,5X .3D.QDQE ,5X .SD.DDOE ,SX ,D . 30E 
NEXT . I 

ASSIGN 3Fat h_4 TO * I CLOSE FILE WITH STANTON NUM6ERS RATIOS 



! CALCULATE LOCAL REYNOLDS NUMBER 

I 

PRINT 

I 

Mu-. 00301 53 ! CINEMATIC VISCOSITY FOR AIR 

Fac=Uinf /Nu 
;< i ' 1 !■ = 1 . 1 5 
X 1 ( 2 ) = i .25 
XI (3)=1 .43 
X 1 \ 4 ) = 1 .50 
XI IS 5-1 .80 
XI <6 >=2.00 



1 880 ! 

1330 FOR 1=1 TO 6 

t 90S Rsy ( I )=X 1 ( I ) *Fac 
1310 NEXT" I 

1320 i 

1 S3© ! CALCULATE AVERAGE STANTON NUMBER 

1940 FOR 1=1 TO S 

1SS0 H‘ I )=S . 

I960 11=1*21-23 

1370 FOR J=M TQ 1+21 

1930 HI ( I >=H! ( I )+H( J ) 

1SS3 NEXT J 

2000 HI < I >=H1 ( I )/2 1 

2313 Still >=H1 < I )/( Ro*Uinf *Cp 5 

2023 NEXT I 

2030 PRINT 

2040 ! 

2350 ASSIGN §Fath7 TO "3TAV2 * 

2360 FOR 1=1 TO 5 

2373 ENTER @Path7 ■, Xt o( I > ,R=o( I ) ,St 1 oil I ) 
2033 Stravv I )=St 1 < I > / 3 1 1 o< I ) 

2030 NEXT I 

2100 CREATE BOAT “STAVFCU" , \ 8 ,8 
2110 ASSIGN §P a t h_5 TO "STnUFCV 
2120 (PRINT AVERAGE STANTON NUMBERS 
2130! 



2U0 PRINT "ROUS X < M > REYNOLDS Nq St/Sto Ave 
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! 073 


i 


1 083 


PRINT "TFLATE - T 1 


! 030 


PRINT 


I 1 33 


PRINT "AVERAGE PLrV 


ms 


PRINT 


1120 


PRINT " NONDIMES I GNi 


I 130 


PRINT 


1 1 40 


PRINT "No X ( M ) 


1 1 £3 


P R I N T ** 


] 153 


PRINT 


1173 


ASSIGN ©Paths TO u : 


1 133 


ASSIGN ©Pat h7 TO 


1 130 


f 


1200 


FOR 1=1 TO 125 


1210 


ENTER ©PathSiHo(I) 


1223 


ENTER ©Fat h 7 ; S t P \ I 


1 233 


NEXT I 


1240 


! CALCULATE THE HEA 


1253 


! 


1253 


r 0 R 1 = 1 TO 


1 273 


Ts I )-T{ 


1253 


: ( : =0/ 


1 233 


St C I ) = H 


1 303 


Str( I )- 


1313 


5 1 f r < I ) 


, 

t jjjU 


NEXT I 


1333 


! ocivcr* rt i c. A rt-AiJ L. 


1343 


2< 1 >=-12.7 


1 333 


X2( 1 >«t . 15 


1 jaw 


FOR 1=2 TO 21 


« *7 T '■*, 

t J i V 


X 2 ( I ) = 1 .15 


1 3S0 


2'. I >=Z(I-i 1+1.27 


13.S0 


NEXT I 


I 400 


ZC 22 >=Z< 1 1 


14-10 


X2( 22 ) = 1 .25 


1420 


FOR 1=23 TO 42 


1 4 J0 


X2( I >=1 .25 


1 440 


Z( I >=Z< I-l >+1 .27 


1 450 


NEXT I 


I 450 


2(43 >=Z( 1 ) 


1470 


X2( 43 >=1.40 


1 4 50 


FOR 1=44 TO 53 


1450 


X2( I >=l .43 


I 500 


Z( I >=Z( I - ! 1+1.27 


1510 


NEXT I 


1 52Q 


Z \ S 4 ) = Z ( 1 ) 


I 530 


*2(64 >=1.5 


« r i n 

l 


FOR 1=65 TO 34 


1553 


X2< I >=i .5 


1 55.0 


Z( I >=Z( I-l) + 1 .27 


1570 


NEXT I 


1 563 


X3= 1 .3 


1 530 


FOR J=5 TO 3 


1 600 


n=J 4 2 1 -23 



z.< cm ) 



H ( U / M * 2 C) St/Sto 



\bf!/ 



t 



17 6 



St /St f " 

* ■* * * 4 ♦ ♦ 






:i50 
2 1 GO 
2 1 70 
2 1 S3 
25 90 
2200 
2210 
2220 
2230 
2240 
2250 
22G3 



2290 

2290 

2300 

2313 



4 o40 
2350 
23S0 



2330 

2333 

2400 

2410 

2423 

2430 

2443 

2453 



2450 

2473 

243G 

2490 

25BS 

2510 

2520 

2530 

2540 

2550 



2570 

2530 

25S3 

2530 

2S10 

2520 

2530 

2540 

2550 

2560 

2570 

2530 



26 S3 



PRINT " * + * + *•**#*■♦*****-**'«■**##****■**-****♦*’•“**■**♦-********* 
FOR r=i 70 S 

PRINT USING 2180? I ,X1 ( I ) ,Rsy< I ) ,Strav( I > 

IMAGE DD.4X .D0.0D.4X .D.4QE.7X ,D.4DE 
OUTPUT QPath_5 ; X 1 < I ) ,Ray{ I ) ,Str*v< I ) 

NEXT I 

ASSIGN 3Pat'n_5 TG « 

I 

GOTO 2SS0 

! SUBROUTINE TO CALCULATE DISCHARGE COEFFICIENTS 



ASSIGN SFathS TO “FILDT10" 

ENTER @Fat n3 5 FIw ,Oelpi 

! GET AVERAGE PLENUM TEMPERATURE 

Tpl=3 

FOR 1 = 123 TO 130 
Tpl=Tpi+T( I > 

NEXT I 
Tpi=Tpi/3 

! CONVERT PLENUM TEMP. TO INJECTION TEMP. 

Tinj-1 . 454G3*To i " . S5S 1 BZ 
Tijab==Tinj+273.15 

Uinj-FIw/. 300921463 IFOR 13 JWJ. HOLES 

!Uinj=FIw/. 0004=5 IFOR 7 INJ. HOLES 

Tor ®T i j ab a - ( U i n j " 2 / ( 2 ♦ 1 3G5 ) / 

R.oc=Panb/ ( 7cr*2S7 / 

Inf l.\=Roc-*Uin j ! INJECTION MASS FLUX 

Pop i=Fanfc/ ( T l j ab 5*2S7 ) ITOTAL DENSITY BASED ON To 

U3— ( Da lp i»2 /Rop i )“ .5 

Pmf i.«,=Roq 1 '•U3 

Cd“Imf ix/rmf 1.* 

Red=Ui nj *. 039525/ . 0G30 1 55 
M=Infix/(Rc*Uinf > • BLOWING RATIO 

Vrat=Uinj/Uinf f VELOCITY RATIO 

Mfr=Roc*Uin j "2/( Ro*Uinf ‘‘2 ) i MOMENTUM FLUX RATIO 



! PRINT 
PRINT 
PRINT 
I 

FRINT 

FRINT 

PRINT 

PRINT 

FRINT 

PRINT 

PRINT 



RESULT S 

“INJECTION TEMPERATURE ( C )=“ i Tin j , "FLENUM TEMP <C> = " » Tp I 
“PLENUM PRESSURE DIFF, « N/M"2 >=” sDslpi 

“COOLANT DENSITY KG/P|*3> = u ?Roc,“ DENSITY RATIO =“ ;Roc/Ro 
"INJECTION VELOCITY < M/S ) = “? Uinj , "MASS FLUX (KG/fTZ SJ-'sImflx 
“REYNOLDS No (Din ) = “ ;Red , “DISCHARGE COEFF , Cd =“;Cd 



FRINT 



FRINT "BLOWING RATIO =“?M ( " VELOCITY RATIO »“ ;Vrat 
FRINT 

PRINT "MOMENTUM FLUX RATIO - "?MPr 

FRINT 

GOTO 590 

MASS STORAGE IS “ : INTERNAL ,4" 

ENu 
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10 

Z0 

3G 

40' 

50 

SO 

70 

30 

50 

I 00 

1 10 

4 ^ a 

t ±VJ 

1 30 
I 40 
t 50 
150 
1 7G 

t n >•% 
i C 

1 33 
230 
213 



REM PROGRAM ACQTrRG 

REM THIS PROGRAM ACQUIRES TEMPERATURES FROM A THERMOCOUPLE 
REM MONTEO ON A AUTOMATIC TRAVERSE MECHAN I SIM AND FROM A FREE 
' ! STREAM VELOCITY THERMCQPLE TO PLOT THE TEMPERATURES PROFILES 
! USING PROGRAM PLQTDE 



! BY ALFREDO ORTIZ , AUGUST ,1987 
REM 



REM VARIABLE NAMES 
REM 



REM E< I ) 
REM- -T( I ) 
REM 



IS THE VOLTAGE READ FROM THE DATA ACQUI3TIGN SYSTEM 
IS THE CONVERSION FROM VOLTAGE TO TEMPERATURES 



REM = AMB=AMS IENT PRESSURE 

REM 7FS=FR£E STREAM TEMP'I MELVIN/ 

REM RG*CENSITY (K3/M'’iJ) 

REM U INF=VbLOC I TV OF THE FREE5TSEAM 
REM 



DIM Y1 i 330) ,Z( 1300 ) ,T! 1000 ) ,T) ( 1000 ) ,D( 1000 > 
REM 



-v — ‘ rs 
U'-'O 


INPUT 


ENTER 


CATE ,(MNGOYY>" ,NS 


Z30 


INPUT 


"ENTER 


T I ME , < HHfiM 1 " ,N3 


240 


REM 






250 


INPUT 


"ENTER 


DISTANCE FROM THE 3.L. TRIP. , 


250 


INPUT 


"ENTER 


POINTS SPANWI3E" ,N.3 


270 


INPUT 


w E N T c. .R 


POINTS VERTICAL (MUST 3E AN IN 


230 


REM N3 


MUST 


oE AN EVEN INTEGER 


230 


INPUT 


"ENTER 


SPANUISE RESOLUTION! IN)" ,24 


300 


INPUT 


u ENTER 


VERTICAL RESOLUTION! IN)" ,Y4 


310 


INPUT 


"INITIAL Z( IN )" Z3 


320 


INPUT 


" INITI 


AL Y( IN ) “ ,Y3 


330 


REM 







543 


Z < 1 )=Z3 


■7f» a 

javj 


Y ( 1 )=Y3 


3S0 


N7=M3/2 


370 


FOR 17=1 


333 


15=17-1 


330 


J 1 = 2 *- 1 5 * 



400 J2 =2 * 1 5 »M3tM3 

413 FOR K- J 1 TO J2 
42G Z ( K )—Z( K- 1 ) + Z4 
43B Y( K ) - Y ( K - 1 ) 

440 NEXT K 

450 J3-2* IS *M3+M3+1 

4S0 Z { J3 )=Z ( J3- i ) 

470 Y( J3 3-Yi J3-1 >rY4 

4S0 J4=J3+1 

493 J5=2»IS*-M3r2*M3 



500 FOR K=J4 TO J5 
510 Z\ K ) = Z( K-1 ) - Z 4 

520 Y < K ) = Y (K-1 ) 

530 NEXT K 
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540 
553 
560 
570 
530 
533 
600 
610 
GZ0 
530 
540 
650 
550 
573 
580 
330 
700 
7 t0 
720 
730 
740 
750 
750 
770 
780 
730 
S00 
810 
320 
S38 
840 
350 
8S3 
a .'0 
380 
3SG 
900 
S18 
920 
330 
940 
950 
950 
970 
933 
930 
1000 
1010 
1020 
1030 
1040 
1 050 
1 050 
1070 



IF I7-N7 THEN 593 
JS-J5+1 
Z(JS)-Z( JG-1 ) 

Y( JS )=Y( J6- 1 >+Y4 

NEXT 17 

REN 

i ENTER FREE 5 7 RE AN TEMPERATURE FRQN DAS. 

OUTPUT 709 s " AI " ; 151 | "UTt " 

ENTER 703 ; X 
E=X* 1030000 

Tf s=- . 0733 1 7 1 - t .026 1 0S45A<^£~ . 000000673 *E*E 
T f s = T f s+273. 1 5 

REN ENTER AMBIENT CONDITIONS 
INPUT "ENTER PANE (IN. OF HG5" ,Panb 

INPUT "ENTER PRESSURE DIFFERENCE (IN GF HZO / " .Del tap 



REM CONVERSION TO 31 UNITS 
Par<ib=Pa.' v )b*33oS . 32 
Rc=Pamb / ( 2S7 »T ? 2 ) 



0eiiap=Deltap*248.7 
REM FREE3TREAM VELOCITY 



r De I ' 



. _ t rv „ 1 ■ 

c3 jj / r.C* ) 



Uinf=( 

REM 

REM ENTER IKE LOOP FOR ACQUIRING EACH TEMPERATURE COMPUTING CGEFF 
REM AND COMPUTING TWO TEMPERATURE COEFFICIENTS AND TUG TEMPERATUR 
PRINTER IS 701 

PRINT CHR3 ( 27 ) ; CHR3 1 38 ) ; CHR3( 1 OG ) ; CHRS ( S3 ) ; CHRS ( 52 ) ; CHRS ( 73 ) 

PRINT CHR3< 27 ) ; CHR3( 33 > ; CHR3( 1 09 > ;CHR5( 49 ) ;CHR3( 78 ) 

l 



PRINT "TEMPERATURE PROFILE COMPUTATION “ 

PRINT 

PRINT "DATE OF RUN IS", NS," TIME =“ ;N3 
PRINT 

PRINT "DISTANCE FROM THE BOUNDARY LAYER TRIP, X 
PRINT 

PRINT "0ENSITY(KG/M-‘3 )“ ,Ro 



PRINT “ FREESTREnM VELOCITY < M/S >* ,Uinf 
PRINT 

PRINT “ PANS I ENT ( N/ M A 2 > " ,Panb 
PRINT 

PRINT "FREE STREAM TEMPERATURE (DEG C) =“ ; Tf s-273 . 1 5 
PRINT 

PRINT "POINTS SPANWISE'' ,M3 
PRINT 

PRINT “F0IN73 VERTICAL" ,N3 
PRINT 

PRINT “SPANWISE RESOLUTION <IN)",Z4 
PRINT 

PRINT "VERTICAL RESOLUTION ( IN)",Y4 
PRINT 

PRINT "INITIAL 2(IN)“,Z3 ; 

PRINT 
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m •-< 



10S0 


PRINT "INITIAL Y< IN) 


iv w *7 
, r 0 




1030 


FRINT 






1 100 


PRINT "* + ** 4 ******-** 




**■**•****♦ + ***+*♦♦**********■♦“ 


1110 


FRINT 






1 120 


PRINT Y 2 




T TF3 DT" 


1 130 


PRINT 


= SSS' 




1140 


FRINT 






1 150 


OISP "HIT <( CONTINUE 


FOR 


DAS )> “ 


1 ISO 


PAUSE 






1 170 


K3=0 






1130 


K2=M3*N3 






1 130 


FOR K-1 TO K2 






1200 


K9-K 






1210 


REM 






1220 


WAIT 10 







1230 REM ACQUIRE THE TEMPERATURE FOR EACH POSITION 
1240 61 =0 

1.253 REM 

12G3 FOR J=1 TO 25 ! ENTER THE DAS ANQ SAMPLE OF TEMPERATURE 25 TIMES 

1273 OUTPUT 733 ; “ AI " ; 1 52 » "VT1 “ 

1230 ENTER 709 sX 
1290 G 1 =0 1 +X 
1 302 NEXT i 
1310 REM 

1320 X-S1/2S ! AVERAGE THE VALUES 

1330 E-X* 1000000 



1 -048 T < K i 335703+ . GZ55E7 1 3*E- . 080300488 *E*fc 
PERATURE 



! CONVERSION FROM VOLTAGE TO 



1353 REM 

1 3SG i GET FREE STREAM TEMP 5 TIMES AND AVERAGE IT 
1 j70 G4=S 



1330 FOR 1 = 1 TO 5 

1390 OUTPUT 733 ; " AI " ; 1 5 1 ; " VT1 

1400 ENTER 703; X 

1410 G4=G4rX 

1420 NEXT I 

1 430 X=S4/5 



1440 E=X * 1 380000 

1450 TICK '- = -.07331 7: r . 02S 1 0S4S4-E- . 008008S78 *E*E 
1 4 o 8 DC K )=T ( K > — T 1 ( K / 

1470 PRINT USING 1 480; Y( K ) ,ZC K ) ,Tv K > , T1 < K ) ,D( K > 

1430 IMAGE MOD . DO ,2X .MOC.OD ,2X ,MDO.DD ,3X .M00.DD.3X , SO. ODE 
1430 REM 
1500 NEXT K 
1510 REM 

1520 MASS STORAGE IS “ : INTERNAL ,4 , 1” 

1530 CREATE BOAT “TPRQ1 4" ,2400 ,3 
1540 ASSIGN OPathZ TO “TPROU" 

1550 FOR 1=1 TO K9 

15S0 OUTPUT 8Path2;Y( I ) ,ZC I ) ,DC I / 

1570 NEXT I 

1530 ASSIGN @Path2 TO 4 

1530 MASS STORAGE IS * ■ INTERNAL ,4 " 

1300 END 



TEN 
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12 i PROGRAM PLOTPRUN 

20 ! THIS PROGRAM READS THE FILE TPRO (AND PLOT TEMPERATURES 

30 ! FROFILES 

40 DIM Y 1 v 330 ) ,Z3( S00 > ,)<£< 300 >C! ] ,01(300 ) ,21 ( 40 > ,22(40 5 ,CS( 13)11] 

50 ! V< I )=D( I ) (TEMPERATURE DIFFERENCE. 

S3 ! 

73 INPUT "ENTER RUN NUMBER " .Runno 
50 FS="DATE=" 

90 MASS STORAGE IS INTERNAL ,4 , 1“ 

1 00 ASSIGN ©Path! TO “TPR07 “ 

110 FOR 1=1 TO S00 

120 ENTER ©Path I ? Y 1 ( I > ,Z0( I ) ,V 1 < I ) 

130 VI ( I >=Yt ( I >*2.54 

140 Z0( I )=Z0( I >'2.54 

150 NEXT I 

150 FOR 1=1 TO 303 



170 


IF 


VI 


( 


T 

X 


K ” . 2 


THEN 


X S \ 


r 

1 


)="0" 








i S3 


IF 


VI 


s 


I 


) \ - . 5 


AMD 1 


VI c I 


) 


>.2 T! 


HEN XT 


A I ) 


« 'M " 

i 


1 30 


IF 


VI 


( 


I 


K = 1 .3 


AMD 


Vi < 


l 


)> . 5 


THEM X3( I 


5 43 •• 2 H 


200 


IF 


VI 


( 


r 


} \ =1,5 


AMD 


VI \ 


I 


)>1 .0 


THEN 


X 3 ( 


I ) * g 3 " 


2 1 G 


IF 


Vi 


( 

A 


I 


>«.© 


AMD 


v 1 { 


I 


) > 1 .5 


THEN 


X3< 


I ) = u 4 u 


220 


IF 


Vi 


( 


I 


K = 2 . 5 


AND 


V i ( 


I 


v r* 

i / L. . 


THEN 


X 5 ( 


I ) = ” 5 " 


233 


IF 


Vi 


f 


I 


) <-=3 . 3 


AMD 


V I V 


I >>2.5 


THEM 


X 5 ( 


I ) - ” £ " 


240 


IF 


VI 


1 


I 


;->3.5 


AND 


V 1 * 


I 


) > 3 . 0 


THEN 


X5( 


I 5= “7" 


253 


IF 


VI 


/ 


V 

l 


;<-4 . 0 


AND 


VI ( 


i 


) >3.5 


THEN 


X 3 ( 


I ; = " 3 " 


253 


IF 


vi 


( 


I 


/ < = 4.5 


AND 


VI ( 


I 


) > 4 . 3 


THEN 


X 3 ( 


I )= " 5 " 


270 


IF 


V i 


\ 


T 

1 


><*5.0 


AND 


VI ( 


r 


) > 4 . 5 


THEN 


X 3 i 


I ) = " d " 


233 


IF 


VI 


\ 


T 


>•'.=5.5 


AMD 


V) < 


i 


) > 5 . 0 


THEN 


X 3 •; 


I )= ,, b M 


230 


IF 


VI 


( 




»S.S 


! nc.'N 


X$( 


7 


)= 11 c '* 









330 NEXT I 

310 6INI7 

320 PLOTTER IS 735, ,, HFGL" 

330 ! PLOTTER IS CRT INTERNAL" 

340 GRAPHICS ON 

353 C3IZE 2.S..S5 

3S0 MOUE 35 , 1 4 

37G LABEL "EMBEDDED VORTEX 10 M/S WITH FILM COOLING, H=0.98" 
380 MOVE 35 , 1 1 

330 LABEL "X = l .480. M , UNKEATEC PLATE .VORTEX GEN. AT 4.75 cm" 
400 MOVE 75,33 

413 LABEL USING "SA . 4 ,60 . 40" ;F5 .Runno 
420 C3IZE 4. 5, -BE 

4 MOVE 43 ,85 

440 LABEL "TEMPERATURE PROFILES" 

453 CSIZE 3.5 , .65 

460 MOVE S7 , 1 7 

470 LABEL "Z (CM) * 

43G DEG 

433 LDIR SO 

530 MOVE 21 ,45 

Si 3 LABEL " Y (CM) " 

520. LDIR. 0 
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530 


3 IEWPQRT 33 ,111 


540 


WINDOW -IS ,5 ,3 


550 


FRAME 


563 


AXES ,=,.5,3,3 


57G 


AXES . 5 , . 5 , - 1 S 


5S3 


AXES .= ,.5,6,1: 


530 


CSIZE 2.0. .65 


600 


MOVE -15.2,11. 


6 i S 


LA3EL " TEMP . D 


620 


A= 1 0 . 65 


630 


FOR K-0 TO 12 


64G 


B“A-K* . 35 


S50 


MODE -14.3,6 


563 


IF K < 2 THEN 


670 


zz< 0 >=-■; z 


633. 


Z 1 < 3 ) = . 2 


533 


22 ( 1 )* . 2 


700 


z i c i i=.s 


7 \ O 


c u, z nz 


720 


Z 2 C K > *? Z 1 C K - ! ) 


73Q 


Z 1 ( K >=Z2C K ) - . 5 


743 


END :r 


750 


IF K/3 THEN 50' 


750 


LABEL USING “D 


770 


GOTO 620 


733 


CSC 10>*‘ a" 


730 


CSC 1 1 > = " b " 


O r* 


CSC 12 )="c" 


8 1 0 


LABEL USING "A 


320 


NEXT K 


— -rs 
dw/o 


CSIZE 1 .3, .72 


340 


FOR 1=1 TO 333 


S53 


MOVE Z 0 < I ) , V 1 < 


350 


IF Y1 C I >>6 .3 Al 


57S 


LABEL USING ”A 


383 


N EXT i 


383 


CLIP OFF 


333 


CSIZE 2.G, .65 


313 


FOR I — .2 TO 1 


923 


MOVE -17.5,1 


930 


I 1 -I + . 2 


S43 


LABEL USING 


3=3 


NEXT I 


363 


FOR J — 17.0 TO 


370 


MOVE J,-.7 


SS3 


J1-J+1 .3 


933 


LABEL USING 


1 3G3 


NEXT J 


1313 


MASS STORAGE I 


1323 


END 
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10 

20 

30 

40 

50 

50 

73 

80 

90 

100 

I 10 

120 

130 

1 40 

150 

150 

170 

180 

190 

200 

213 

220 

230 

240 

250 

2S0 

270 

230 

290 

300 

310 

320 

330 

540 

350 

350 

370 

380 

390 

400 

410 

420 

430 

443 

450 

453 

470 

480 

490 

500 

510 

520 



! PROGRAM PT3TAV 

(THIS PROGRAM READS A DATA FILE AND PLOTS AVERAGE STANTON NUMBERS 
I VERSUS REYNOLDS NUMEERS. 

1 

DIM Rey<6 ) .X 1 <5 ) ,Stl (6 ) .CSC SI .St <6 5 
MASS STORAGE IS *: INTERNAL ,4 , 1 " 

ASSIGN §Pat h2 JO "STAV1" 

FOR 1=1 TO 5 

ENTER SPathZiXH I ) ,Rey< I ) ,St 1 < I ) 

NEXT I 

PRINTER IS 1 

PRINT “ENTER RUN NUMBER TO WHICH THIS DATA CORRESPOND MMDDY Y . HHMM > “ 
INPUT Runno 
F$ = " DATE = " 

GINIT 

PLOTTER IS 70S,“HPGL" 

• PLOTTER IS CRT INTERNAL” 

GRAPHICS ON 
CS I ZE 2.5 , . 55 
MOVE 42 , 14 

LABEL “FEEESTREAM IS M/S, NO FILM COOLING - 
CSIZE 4. 5. .55 
MOVE 40 ,95 

LABEL “AVERAGED STANTON NUMBERS “ 

CSIZE 2.5..SS 
MOVE 75,53 

LABEL USING “BA ,t ,S0 . 40“ ; F$ , Runno 

CSIZE 3. 4,. 55 

MOVE 52,17 

LABEL "RE* * 10 "6" 

DEG 

LOIR 90 
MOVE 22 ,33 

LABEL "STANTON Mo *10" -3" 

LOIR 0 

VIEWPORT 30,1 1 S ,25 ,S3 
WINDOW 0,5, 0,5 
FRAME 

AXES .5 , .5 ,0 ,0 ,2 , 1 
AXES .5,. 5, 5,5,2, 1 
CLIP OFF 
CSIZE 2.0 , .65 
LORG 2 

FOR 1=0 TO 6 STEP .5 
MOVE -.32,1 

LABEL USING " # ,D.DO“ ; I 
NEXT I 
LORG 4 

FOR J =0 . TO 5.0 STEP .5 
MOVE J ,-.22 

LABEL USING “#,DD.D";J 
NEXT J 
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530 


! ENTER VALUES FOR ANALYTICAL FORMULA 


540 


Z l = 1 . 1 i METERS 






550 


Pr= . 7 1 IPRANOTL 


No FOR 


AIR 


550 ' 


FOR 1=1 TO 5 






570 


£'t( I ) = < ( .030*( Rey( I ) )'•( - 


.2 ) >*< 1 , 


Zi/XI < I 


530 


S.t ( I )=St(T)*l .0E+3 






530 


Rsy.( I )=Ray < I )*\ .0E-S 






600 


PLOT Rev < I ) ,St( I ) 






SI0 


NEXT I 






520 


CSIZE 1.3, .72 






630 


FOR 1=1 TO S 






540 


CSCI 3 = ''a" 






550 


St!( I >=St 1 ( I ) * 1 . 0Ei-3 






650 


MOVE F-ey ( I ) ,St1< I ) 






570 


LABEL US IMG "A''iC5CI] 






530 


PENUP 






530 


NEXT I 






700 


HASS STORAGE IS 11 : INTERNAL ,4 M 




710 


ENO 
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10 ! PROGRAM PTSTLC 

20 'THIS PROGRAM READS A DATA FILE AMO PLOTS SPANWISE VARIATIONS 
30 ! OF HEAT TRANS FER COEFFICIENTS 

40 ! 

50 DIM H< 200 > , St < 230 ) ,Rey( 5 ) ,2( 200 > ,X 1 ( 6 ) ,X2< 200 > ,B( 1 26 > ,C$< S >C 8 

60 MASS STORAGE IS INTERNAL ,4 , 1" 

70 ASSIGN @Pat h 1 TO “HDAT2" 

80 FOR 1=1 TO 126 

90 ENTER gPathi ;H( I > ,St < I > ,X2(I > ,Z( I > 

100 NEXT I 
110 F$="OATE =" 

120 INPUT -ENTER RUN NUMBER ( MMOOYY . HHMM ) " ,Runno 
130 6INIT 

140 PLOTTER IS 705 , '‘HPGL" 

150 ! PLOTTER IS CRT INTERNAL" 

1S0 GRAPHICS ON 

170 CSIZE 2.5 ,.65 

1 30 MOVE 42,14 

130 LABEL "FREE STREAM 10 M/S, NO FILM COOLING" 

2O0 CSIZE 4 . 5 . . 65 

210 MOVE 20,55 

220 LABEL "SPANUISE HEAT TRANSFER COEFFICIENT" 

230 CSIZE 2 . 5 . . SE 

240 MOVE 75,33 

250 LABEL USING “6A ,* ,SD.4D" ;F$ .Runno 

2S2 CSIZE 3. 5, .65 

270 MOVE 55,17 

280 LABEL "Z ( CM ) “ 

290 MOVE 22 ,35 

300 DEG 

310 LOIR 90 

320 LABEL "H U/M A 2 DEG C" 

330 LOIR 0 

340 VIEWPORT 30 , 1 IS ,25 ,83 
350 WINDOW -15.24,15.24,0,130 
360 FRAME 

370 AXES 1 .27,10,-15.24,0,2,2 

380 AXES 1.27,10,15.24,100,2,2 

390 CSIZE 2.0,. 55 

400 MOVE -13,23 

410 LABEL “X < M ) " 

420 A=22 

430 FOR K-l TO 6 

440 0=A-K«3 

450 C$( 1 5="=X=I . 15" 

460 C$<2 )-"*X-l .25" 

470 CSC 3 >="=X=! .40“ 

480 C$( 4 )="=X=1 .60" 

490 CS(3)="=X=I .80" 

500 C3( 6 ) = - =X=2 . 00 " 

510 MOVE -1 4 ,0 

520 LABEL USING “0 ,8A“ ;K ,C5? K >C 1 ; 3 3 
530 PENUP 

540 NEXT K 
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553 CSIZE I .3, .72 
550 FOR J-l TO 5 
570 M=J *2 1 -20 
5S0 FOR I-M..TO J*2 I 
5S0 B( I >=J 
600 MOVE Z < I ). , H ( I > 

610 LABEL. .USINS “D";B(I) 

520 PENUP- 
630 NEXT I 

643 NEXT J' 

653 CLIP OFF 

653 CSIZE 2. 3, .55 
570 LORG 2 

580 FOR 1=0 TO 103 STEP 10 
SS0 MOVE -17.75,1 

703 LABEL USING "* .DDO.D" i I 

•710 NEXT I 

720 LORG 5 

730 FOR J— 15.24 TO 15.2* : _ EP 5 
740 MOVE J ,-2.0 

753 LABEL USING “4 .ODD .DD" ; J 

750 NEXT J 

770 /MASS STORAGE IS INTERNAL ,4 
780 END ’ 



.08 



186 



10 


•PROGRAM FL3TRTI0 




20 


•THIS PROGRAM READS A DATA FILE AND PLOTS 


THE RATIO OF 


33 


• DOER BASE LINE STANTON NUMBERS VERSUS REV 


HOLDS NUMBERS 


40 


j ...... 




53 


DIM Re y ( S > , X i ( G > , 5 1 H S > , 5 1 Z ( 5 ) , Re 1 ( B ) , C$ i S 


>[ 1 1 ,PS(B)C 1 1 


50 


DIM S 1 3 ( B ) ,fte 3( S ) ,St 4< S ) ,Re4CS ) ,Strl IS > ,Str2(S > ,Str3< B > 


70 


ASSIGN ®Fath1 TO "STAv'" 




30 


FOR 1*1 TO E 




90 


ENTER SPathl 5 XUI >,ney< I 1 ,Si1< I > 




130 


NEXT I 




1 13 


ASSIGN ®Path2 TO “STAUFC" 




120 


FGR 1= ! TO 5 




; 30 


ENTER ®Path2 ? Re i ( I ) ,Si2( I ) 




149 

i C ."‘i 


str ! 1 1 > = s t z i ;/st mi; 

MU' * T 




1 SU 

160 


ricA I i 

ASSIGN §Pst h3 TC "STAUF1 




! 7G 


FOR 1=1 TO S 




! 33 


ENTER SPatn3«Re3i I ' ,3t3< I > 




1 S3 


3tr2 \ I ) = 5 1 3 ( I ) / S 1 1 \ I J 




233 


NEXT I 




213 


ASSIGN | 5Path4 TO "STAUF2" 




223 


FOR 1=1 TO 5 




230 


ENTER § P ath4 ;Re4i I ) ,SU( I > 




240 


S t r 3 1 1 > = 3t4‘. I >/ St! ( I 5 




L 3 4- 


NEXT I 




253 


PRINTER 13 1 




~>"v. 

U i 


PRINT “ENTER TODAY’S DATE ( MMODYY .HHMM ; ' 




232 


INPUT Run no 




230 


F 3= "DATE = " 




330 


PRINT “ENTER IN ORDER CORRESPONDING BLOWING PATIOS, Ml, 


3 ! 9 


FOP. 1 = 1 TC 3 




—fr^rx 


INPUT 11 ( I ) 




330 


1 ECHO PRINT M 




340 


PRINT "M=" ;M( I ) 




353 


SINIT 




350 


! PLOTTER IS 70S,"HPGl" 




370 


PLOTTER IS CRT, "INTERNAL." 




3S0 


GRAPHICS ON- 




390 


X_gdu_fia /:= r CO ‘MAXt 1 , RATIO) 




433 


Y_gdu_i'ia.-: = 1 00+MAX t 1 , 1 /RATIO 1 




410 


CSIZE 4.5, .55 




4Z0 


LORS S 




433 


MOVE X_gdu_na.*:/2 , . 9-»Y_cdu_r , ia>: 




440 


LABEL "FILM COOLING* 




450 


DEG 




463 


LOIR 30 . 




470 


CSIZE 3.5 




433 


MGUE . 1 ♦ X_gdu_ma.-. ,Y_gdu_ma:-i/2 




493 


LABEL "St /St 0 * 




500 


LOIR 0 




510 


LOR 6 4 




523 


MGUE X_gdu_:na.-; /2 , . 20-»Y_gdu_ina.<. 




533 


LABEL "REYNOLDS No* !0''S" 
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543 CSIZE 2.5 
350 LGRG 3 

5S0 MOLE . 3*X^du,jnax , . o5*Y_gcfti_tta>; 

570 LABEL USING "BA ,4 ,SD . 40° ;F3 ,Runno 

583 VIEWPORT . Z *X_gdu_Max 9 *X_gdu_>ax , . 30*Y_gdu_max , . 30*Y_gdu_p3 
590 FRAME 

500 WINDOW .5,1 .5 ,.3, 1 .0 
5 » 0 AX E 3 . 1 , , 1 , . 6 , . 3 , 1 , 1 

S23 AXES .1 , J ,1 ,S ,1 .3,1 ,1 

530 CSIZE 2.5 
540 MOVE 1 .3 ,.5 
630 LABEL "BLOWING RATIOS " 

55 3 A= , 55 

570 FOR k-j TQ 3 

550 u=A-K * . 2 

533 33 ( 1 5= " ***«*■’ 

700 59^2 ' - ’’ •;> =N- " 

7 ■ 3 33(3;-“ " 

720 MOLE i .28,0 

73S LABEL USING "4A ,$ ,0.00" »B$< I) ,M( I ) 

7-0 NEXT K 

733 CLIP OFF 

TBS CSIZE 2.3.. 55 
770 LOR 5 2 

733 FOR I -.3 TO 1 .0 STEP . 1 
750 MOLE .54,1 

SGO LABEL USING ‘'£,D.D M *I 

313 NEXT I 

32Q LORG S 

530 FOR J = ,G TO I *.7 STEP . t 

340 MOvE J , . 27 

S50 LABEL USING "*,Q.O" :J 

SEC NEXT J 

873 CSIZE » . 3 , » 72 

830 FOR 1 = 1 TO 6 

S30 CSC I 

900 MGUE ReyC I > .Sir I ( I > 

310 LABEL USING " A"iC3(I> 

333 S3 ( I ) = " o ,r 

933 MOVE Reyv X > ,Str2C I > 

940 LABEL USING "n";G3(I) 

55G P$< I >="+" 

363 MOUE ReyC I ) ,Str3( I ) 

373 PENUP 

353 NEXT i 

330 END 
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10 ! PROGRAM PLSTRVOR 

20 'THIS PROGRAM READS A OATA FILE AND PLOTS SPANWISE VARIATIONS 

30 ! OF STANTON NUMBERS RATIOS BY ROUS 

40 ! 

50 DIM St r (200) ,Rey( 5 ) , Z( 200 ) ,XI!8 ) ,X2( 200 > ,C'S< S > C 3 1 ,B< 1 2B ) 

60 MASS STORAGE IS INTERNAL ,4 , 1" 

70 ASSIGN ©Path! TO “STRDT2 " 

S0 FOR 1-1 TO 126 

30 ENTER ®PathliStr(I),X2<I),Z(I) 

100 NEXT I 

110 F£=‘‘DATE =" 

120 PRINTER IS 1 

130 PRINT "ENTER RUM NUMBER TO WHICH THIS PLOT CORRESPOND ( MMDDYY . HHMM ) 
1 40 INPUT Ruinno 

150 GINIT 

160 PLOTTER IS 70S ."HPCL" 

170 ! PLOTTER IS CRT INTERNAL" 

1S0 GRAPHICS ON 

I S3 CSIZE 2.5, .55 

200 MOVE 42,14 

210 LABEL "FREE STREAM 10 M/S, NO FILM COOLING " 

220 CSIZE 4. 5,. 65 

230 MOVE 23, 85 

240 LA8EL "STANTON NUMBER RATIOS, VORTEX GEN 42" 

250 CSIZE 2.5 , .55 

2E0 MOVE 70,33 

270 LABEL USING "SA ,4 ,60 . 4D” «F3 .Runno 

2S0 CSIZE 3. 5, .65 

230 MOVE 65,17 

300 LABEL “Z (CM) ‘ 

310 MOVE 22,45 

320 DEG 

330 LDIR 30 

340 LABEL “St/Sto" 

350 LOIR 0 

3S0 VIEWPORT 30,116,25,83 

370 WINDOW 23, -20,. 4,1 .4 

3S0 FRAME 

330 AXES 5 , . 1 ,20, .43,1 ,1 

400 AXES 5 , . 1 ,-20 ,1 .4,1,1 

410 CSIZE 2.0, .55 

420 MOVE 13. 5,. 85 

430 LABEL "X (M) u 

440 A= . 80 

450 FOR K= I TO 6 

460 D=A-K* . 04 

470 C$( 1 )="=X=1 . 15” 

480 C$(2 >="=X=1 .25" 

490 C3< 3 >="=X=1 .40“ 

500 C3(4)="=X = 1 .60" 

510 C$<5/-"*X-1 .80" 

520 C$( 6 ) = ” =X=2 . 00” 

530 MOVE 14,0 
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540 

5 = 0 

560 

570 

580 

590 

600 

510 

S20 

630 

640 

550 

680 

670 

680 

693 

700 

710 

720 

730 

740 

753 

750 

773 

780 

790 

80S 



LABEL USING “D ,8A" ;K ,CS< K ) 

PENUP 

NEXT K 

CSIZE 1.3,. 72 
FOR J=1 TO 6 
M= J+2 1 -20 
FOR I-M TO J *2 1 
B( I )=J 

MOUE Z( I ) ,Str< I ) 

LABEL USING "D” ;B( I ) 

PENUP 
NEXT I 
NEXT j 
CLIP OFF 
CSIZE 2.3, .55 
LCRS 2 

FOR 1= . 4 TO 1.4 STEP . 1 
MOUE 21.7,1 
LABEL USING "4,0. O' 1 : 1 
NEXT I 
LORG 6 

FOR J =-20 TO 20 STEP 5 
AGUE J , . 37 

LABEL USING " Z ,DQD .DO" ; J 
NEXT J 

MASS STORAGE IS INTERNAL ,4 " 
END 
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to 

20 

30 

40 

50 

60 

70 

90 

93 



! PR06RAM PL3TRUU6 

1 THIS PROGRAM REAOS A DATA FILE ANO PLOTS SPANWISE VARIATIONS 
! OF STANTON NUMBERS RATIOS BY ROWS FOR FILM COOLING DATA ANO 
IFOR VORTEX OATA 

t 

DIM St r 1 < 1 2S > ,Str2( ! 2S ) ,Z< 1 25 ) ,X2 (i ZS > ,X3( 1 2S > ,C«( 1 26 >C 1 3 ,B5( { 2S >C 1 1 
DIM St fi 1 2b >,ZK 12B) 

| 

MASS STORAGE IS * : INTERNAL ,4,1* 



130 


ASSIGN gPathl T 


1 5 3 


ASSIGN ? 3Patn2 I 1 


120 


f 


1 3® 


FCF, 1 = 1 TO 125 


1 43 


ENTER S s 3 1 h ] ; 5 1 : 


I 5G 


ENTEP OP a t h 2 : 5 u 


! BO 


NEXT I 


170 


F 5=“ DATE =" 


1 


PRINTER IS > 


! 20 


PE [NT “ENTER PUi 


202 


INPUT R.innc 


2 ! 0 


C-INIT 


220 


FLOTTEF IS 735, 


230 


'PLOTTER IS CRT 


2 40 


GRAPHICS ON 


250 


C3I2E 2.5 , .65 


2SG 


MOVE 26,14 


270 


LABEL "FREE STR 


Zb'S 


CSI2E -i-.S , .EE 


290 


MOVE 43 ,65 


330 


LASEL “STANTON 


313 


CSIZE 2 .5 , . E5 


320 


MODE 75,93 


330 


LABEL USING “SA 


340 


CEI2E 3.5. .55 


350 


MOVE 55,17 


350 


LABEL '2 (CM 5 


372 


MOVE 25,45 


383 


DEG 


3S3 


LOIR 93 


402 


LABEL ‘ £ t / 5 1 o ‘ 


413 


lDIR 3 


420 


VIEWPORT 36,516 


433 


WINDOW -20,20,. 


440 


FRAME 


450 


AXES 5 , . 1 .-20 , . 


463 


AXES 5 , . i ,20,1. 


473 


CSIZE 3.0. ..65 


480 


MOVE -17.5,1.0 


490 


LABEL "K * 2.00 


530 


CSIZE 2. 5. .55 


510 


MOVE 0. , .4 


52 B 


LABEL "A =3 . 47" 


530 


MODE -17. 5,. 41 



CORRESPOND ( MADE V ? . HUM. A } ' 
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54.0 


LABEL " o = WITHOUT VORTEX 


550 


MOVE -17, 5,. 37 


553 


LABEL “ + = WITH VORTEX" 


573 


CSIZE 1 . 3 , .72 


583 


FOR 1=106 TO 1ZS 


530 


8 $< I ) = " g " 


S 33 


CSv I ' * :l + “ 


513 


MOVE I( I ) ? Str H I ) 


520 


LA 6 E ; _ USING “A" :S3(I) 


520 


MOVE Z T \ I ; .StrZt I ) 


54Q 


LABEL USING “ A ” ; C 3 v I ) 


5 So 


-EfsiUF 


r 0 
c _• v; 


NEXT I 




CLl? Or ^ 


GS0 


CSIZE Z.U..5E 


r* n •“» 
C 7 'u 


lCF.G Z 


720 


FOR I = . 3 f 0 i . 2 STEP . 1 


/ 1 >: 

7Z0 


LABEL USING *’2,0.0" ;I 


730 


NEXT I 


7-^-0 

7C;7 


LORG S 

rpc : n t n 7 .*?. 7 t -x o c 


750 


MOVE J , .27 


770 


LABEL USING ” $ f DOD . DO " 


73S 


N E a T J 


7 SO 


MASS STORAGE IS ” : INTERNA 


533 


END 
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IQ 


! FROGRi 


20 


l THIS 1 


30 


!GF ST 


43 


t 


50 


i ORTIZ 


60 


DIM St 


70 


MASS S 


80 


A5SI5N 


33 


PRINTEi 


103 


FOR 1 = 


113 


ENTER i 


i 20 

4 7 ,n 


NEXT I 


1 ^ <J 

140 


PRINTS! 


153 


PRINT 


153 


INPUT i 


170 


) 


! 53 


5 1 - i r - 


130 


L i . 


r t •* 7 


FOR 1 = 


I 1 V? 


T F 5 x r 


222 


i “ c t 


413 


IF S t r 


41 i 


IF Str- 


412 


IF Str 


4 T C 


NEXT I 


440 


i 


450 


SIN IT 


4 50 


FLO T TEi 


470 


1 PLQTT! 


43G 


GRrtPH I r 


490 


C b i Z b 


500 


MOVE 1. 


510 


LABEL 


520 


CS 1 IE 


530 


MOVE 4 


E40 


LABEL 


550 


CSIZE : 


550 


MOVE 7' 


573 


LABEL i 


C33 
w w W 


CSIZE 


550 


MOVE 5 


500 


wAEEL 


£10 


MOVE i 1 


520 


CES 


523 


LDIR Si 


840 


LABEL 


S52 


LDIR 0 


553 


cbize : 


873 


MOVE 5' 


580 


LABEL 



^ UUMN| ' * 



3Q 



I '.J > b 



. b , . bb 
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SS0 


A=24 




730 


FGR K=0 TO 2 




7 ! 3 


Q-Fs-K *2.5 




720 


MOVE 34,5 




730 


Z2 i 0 )=S tni n-:-K * Inc 




740 


Z ! ( 3 ) = S c i n + < K + ! ) » I r. c 




750 


Z2( ! .'=.38 




750 


ZH 1>=! .02 




773 


:2(2>=1 .32 




730 


Z ! i 2 1 = 1 .10 




750 


LAG EL USING 30Q:K ,Z2< K > , 


ZHK 


see 


IMAGE D , 2 X .0 . ZD ,2X ,0.30 




p 1 


NEXT K 




SZ0 


-OR r.-3 TO 4 




330 


5=fW K - 3 ; ■> 2 .5 




840 


MOVE 54,0 




O !d w* 


Z 2 ( 3 > = 1 . 1 3 




sse 


Z! Z i .20 




370 

O J r ^ 


22 . 4 :• = ! .27 






2. : ’ T ) - ! . _• .v 

LABEL USING 800 : K ,22 ( K ) , 


7i (i.; 


300 


NEXT K 




310 


') I EL? GET 25,1 U ,40,70 




3 Z 0 


U'lMOG'L ! . !3 ,2,20 ,23. ,-20 


, 


cu 

u 

C2 


FRAME 




348 


AXES .35 ,5. , 1 . ' 0 ,23. ,4 , 1 




350 


AXES .05,5. ,2.23,-23. ,4, 


i 


553 


rcCVE i . 33 , i 7 




370 


LABEL " M =3 . SS " 




550 


■~CT->C 1 T ■?-> 

•- w i _ w 1 • - , • • 4- 




530 


FOR 1=1 TO i 26 




1030 


MOVE X 2 ( I > ,2(1 > 




1 0 1 G 


LABEL USING ‘ D * : 5 1 I > 




IflZC 


PENL'r 




1030 


NEXT I 




i 040 


CLIP OFF 




1 050 


CSIZE 2.0,. 55 




1050 


L0R5 2 




1 Q7G 


FOR I =-23 TO 23 STEF S 




! 030 


MOVE 1 .35 , 1 




1030 


LABEL US I NG “ 5 , ODD . 0 “ s 


K 


1 130 
1 ! ; 0 


NEXT I 
LORO- 3 




l 1 20 


r OF J=i . tO TC 2.33 STEF 


. 1 


1 133 


MOVE J ,22 




i 1 40 


LABEL USING “4 ,3.00" ; 


J 


1 150 


NEaT J 




1 ! GG 


MASS STORAGE IS *: INTERNAL ,4 


i 170 


END 





L- 
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10 ! PROGRAM FLSTRFC 

20 ! THIS PROGRAM READS A DATA FILE AND PLOTS SPANUISE VARIATIONS 

30 ! OF STANTON NUMBERS RATIOS BY ROUS FOR FILM C00LIN6 DATA 

40 ! 

50 DIM StrC 200 ) ,Rey< 6 ) ,1(200) ,X 1 < 5 ) ,X2< 208 > ,CS( 6 ) C 8 ] ,B( 12S> 

GO MASS STORAGE IS *: INTERNAL ,4 , 1 * 

73 ASSIGN SPatn! TO “STRFC" 

20 FOR 1-1 TO I2E 

90 ENTER SPathl ;Str< I ) ,X2( I ) ,Z( I ) 

! 30 NEXT I 

1 ! 0 F$='' DATE =” 

! 20 PRINTER IS I 

130 FRINT "ENTER RUN NUMEER TO WHICH THIS FlOT CORRESPOND < MMODYY . HHMM ) 
140 INFL'T Runnc 

150 -3 INI T 

ISC PLOTTER IS 70S , “ HF'Gl “ 

170 I PLOTTER IS CRT , " INTERNAL" 

ISO GRAPHICS ON 
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r* -z T *• T. ”> C 2C 

lrCi.-L :.r ,.33 


1 0 0 


MOVE 4C- ,14 


2 1 ‘0 


LABEL "FREE STR 1 


220 


C3I2E 4. 5 ,.55 


n-rr* 


MOVE 42 .35 


240 


lABcl. 'S ! ANTON 


2 5 Vw 


C'SIZE 2.5,. S3 


^ p ™ 

4. Z, <u 


MOVE 75,33 


:?e 


LABEL USING “SA 


2B<0 


C S I Z E 3.5, . S E 


4— W 


MOVE 55,17 


333 


LAEEL “2 (CM) 


313 


MOVE 25,45 


320 


DEG 


330 


LOIR 30 


340 


LABEL “St -'St o' 


330 


lDIR q 


yr ,> 
w tl *J 


VIEWPORT 30,1 IS 


370 


WINDOW -23,20,. 


330 


FRAME 


332 


AXES 5,.1,-ZO,. 


400 


AXES 5 , . 1 ,20 ,1 . 


410 


CSIZE 2.0..E5 


420 


MOVE -17.5 , . 65 


430 


LABEL "X 1 M 5 “ 


440 


A*. 52 


450 


FOR *■' = ! TO 5 


460 


D=A-K ♦ . 24 


470 


CSC 1 > = “=;< = ! .15" 


4S0 


C S ( 2 ’= “ =*X= 1 .25“ 


493 


C3<3 )="=X-i .40" 


530 


C$< 4 .’*=“=:<«= i.GO" 


5 1 G 


C$( 5 )=" =X=\ .80“ 


320 


CSC 6 )="='/.=2.QQ“ 


530 


MOVE -13. ,D 


540 


LAEEL USING “D, 



TREAM 5 5 M/S, WI- 



LL M COOLING 
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c 



550 


PE MOP 


560 


NEXT K 


570 


CSIZE 1 .3,. 


5 = 0 


FGR J* ! TO 


550 


-20 


SOS 


FOR I =M TO 


612 


E( I )-J 


620 


MGUE Z ( I ) , S 


53G 


LABEL USING 


5^0 


PENUP 


550 


NEXT I 


550 


NEXT J 


872 


CLIP OFF 


S80 


CSIZE 2.3 ,.! 


59G 


LORG 2 


733 


FOR I-.2 TO 


*■“' ♦ rx 
( i W 


MOVE -21 . 


7 n 0 


Lr.wSL LSI 


/ -■ 0 


NEXT I 


7i0 


LORG 5 


753 


FOR i=-20 T 


753 


NOvE J , . 


773 


LASEl US 


790 


NEXT J 


753 


NASS STOP AG 


833 


END 



) 



:E< 



) 



- 



,L ; 



STEP : 



't ,CDD 
,f : I NT 



i 

I 



nr " . 
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APPENDIX D. 



DATA FILES CATALOG 



n C rV 



’T T > . . . . 
i i. 



-'^ER DAT; 



DISk No. 



No. 


File Name 


Dale 


1 iijee. 


v or Lex 


Description 


001 


TDAT 1 


71787.1702 


XT 

X N 


N 


15 m/s, baseline 


002 


ID ATI 


" 


N 


M 


" 


003 


STAV1 


" 


N 


N 


baseline Si av 


004 


HDATl 




N 


N 


baseline local h 


003 


TDAT2 


72087.1315 


N 


N 


10 m/s base line 


006 


IDAT2 


" 


N 


N 




007 


STAV2 




N 


N 


baseline Slav 


00? 


II DAT 2 




N 


NJ 


baseline local ft 


009 


STRDT2 


72187.1524 


N 


y 


10 m/s St/Sto 


010 


STRV2 


1 ■ 


N 


Y 


10 m/s (St/Sto ) a v 


Oil 


TDAT3 


1 ' 


N 


Y 


T(I) 


012 


IDAT3 




N 


Y 


input data 


013 


TDAT4 


72187.1748 


N 


Y 


15 m/s Till 


014 


IDAT4 


'* 


N 


Y 


input data 


015 


STRDT1 


' 1 


N 


Y 


St/Sto 


016 


STARV1 


M 


N 


Y 


(St/Sto) a v 



017 


TDAFC2 


72387.1318 Y 


N 


10 m/s, M - 0.68 T(I) 


018 


IDAFC2 


Y 


N 


input data 


019 


FILDT2 


Y 


N 


film cooling data 


020 


STRFC2 


Y 


N 


St/Sto 


021 


STAVFC2 


Y 


N 


(St/Sto )av 
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022 


ID AFC 1 


72387.1622 


V 

i 


N 


15 m/s. Til) 


023 


I D AFC 1 


" 


V 

1 


\T 
* « 


, input data 


024 


FILDT1 


* 


V 

i 


N 


,lnj. data 


025 


STRFC1 




Y 


N T 


. St/Sto 


026 


STAVFCl 




Y 


N 


, 1 St/Sto i-av 


027 


TDAFC4 


72487.1401 


Y 


N 


10 m/s, repeat 


028 


IDAFC4 


" 


V 

l 


V 

lN 




i . ^ A 


rriT at ; 










v/_ ; 


X 1 I 




'■ 






u j n 


STRFC i 


" 


T » 


*- r 




/A i 

u J i 


STAVFC4 




\~ 

l 


V 




052 


TDAFC3 


72487 1621 


Y 


N 


1 5 m/s. repeat 


053 


IDaFCF 




V 

1 






034 


FILDT3 




Y 


\I 




035 


ST RFC 3 




\r 

I 


M 




036 


STAVFC5 




Y 


N 




037 


TDAFC6 


72587.1401 


Y 


N 


10 m/s. M = 1.26 


038 


IDAFC6 


1 ' 


Y 


\i 




039 


FILDT6 




Y 


N 




040 


STRFC6 




i 


N 


. St/Sto 


041 


STAVFC6 




Y 


N 


. lSt/StO)av 


042 


TDAFC5 


72587.1808 


Y 


N 


15 m/s, M = 0.86 


043 


IDAFC5 




Y 


N 




044 


FILDT5 




Y 


N 




045 


STRFC5 


“ 


V 

1 


N 


.St/Sto 


046 


STAVFC5 




Y 


M 


, 1 St/Sto i av 


047 


TDAFC7 


72887.1206 


Y 


N 


15 m/s. M = 0.86 


048 


IDAFC7 


" 


Y 


N 


, Repeat 


049 


FILDT7 


" 


Y 


N 




050 


STRFC7 


" 


Y 


N 




051 


STAVFC7 




Y 


N 




052 


TDAFC8 


72887.1353 


Y 


N 


10 m/s, M = 0.68 


055 


IDAFC8 




Y 


N 


. Repeal 


054 


FILDT8 


" 


Y 


N 




055 


STRFC8 


" 


Y 


N 


, St/Sto 


056 


STAVFCS 


•* 


Y 


N 


, (St/Sto i a v 
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057 


STAFC8 


Y 


N 


Stf numbers. 


058 


TDAFC10 


72887.17 42 Y 


V 

l 


10 m/s. M - 0.68 


059 


ID AFC 10 


Y 


Y 


Vortex sen. 4.79 cm 


060 


FILDTIO 


Y 


Y 




061 


STRFCV 10 


Y 


Y 


St/Sto 


062 


STFCV 1 0 


Y 


Y 


lbt/Stoi a v 



065 


TDaFCI 1 


72987.1508 


Y 


V 

L 


t c m/s. Ai - u.'t 


064 


ID AFC 11 


1 


Y 


Y 


Vortex gen. <<■•< 4.79 cm 


065 


FILDT1 1 




Y 


V 

i 




0 6 o 


STRFCV 1 1 




Y 


Y 


St/Sto 


067 


STAFCV 1 1 




V 

l 


Y 


i St/Sto J a v 


068 


TDAFCI 3 


72987.1345 


Y 


V 

i 


1 5 m/s, M - 0.47 


069 


ID AFC 1 3 


1 1 


Y 


Y 


Vortex gen. & 3-52 cm 


070 


FILDT13 


" 


Y 


Y 


" 


071 


STRFCV 13 


• ' 


Y 


Y 


' 1 


072 


STAFCV 13 




Y 


Y 




073 


TDAFCI 5 


72987.141 5 


\r 

I 


Y 


15 m/s, \1 = 0.47 


074 


ID AFC 1 5 


“ 


Y 


Y 


Vortex gen. <?> 6.06 cm 


075 


FILDT15 


'* 


Y 


Y 


" 


076 


STRFCV 15 




Y 


Y 


1 1 


077 


STAFCV 15 




Y 


Y 




078 


STAFC6 


72587.1451 


Y 


N 


10 m/s, M' 1 .26. Stf No 


079 


STAFC7 


72887.1206 


Y 


N 


15 m/s. M=0.86, Stf No 


080 


STAFC3 


72487.1621 


Y 


N 


1 5 m/s, M=0 47. Stf No 


081 


TDAFC21 


80187.1354 


Y 


Y 


15 m/s, M = 0.86 


082 


IDAFC21 


*' 


Y 


Y 


Vortex (a 6.06 cm 


083 


FILDT21 


" 


Y 


Y 


” 


084 


STRFCV2 1 


" 


Y 


Y 


St/Sto 


085 


STFCV 21 




Y 


Y 


averaged St/Sto 


086 


TDAFCI 7 


80187.1507 


Y 


Y 


15 m/s, M« 0.86 


087 


IDAFC17 


“ 


Y 


Y 


Vortex & 4.79 cm 


088 


FILDT17 


'* 


Y 


Y 


• i 


089 


STRFCV 17 




Y 


Y 


" 


090 


STFC17 


" 


Y 


Y 
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09 i 


TDAFC19 


80187.1440 


Y 


Y 


15 m/s, M= 0.86 


0°2 


IDAFCI9 


" 


Y 


V 

i 


Vertex 3.52 cm 


093 


PILOT 19 


' 


Y 


Y 


" 


09 « 


ST RFC V 1 9 




V 

l 


Y 




09 3 


STFCV19 




Y 


Y 




096 


TDAFC12 


80187.1841 


Y 


T T 
1 


10 m/s. M = 1.26 


097 


ID AFC 12 


" 


Y 


Y 


Vortex (S> 4.79 cm 


, a n o 

\J V o 


PILOT 12 




v 

1 


1 




non 


Z> i i\iK. V I - 


“ 


* r 

i 


T T 

i 


St/Sto 


100 


STFCV 1 2 




AT 

1 


\r 

i 


Si/S‘c 


101 


TDAFC14 


80487 1321 


Y 


Y 


10 m/s, AW1.26 


102 


IDAFC14 


'' 


Y 


Y 


Vertex (?) 4.7° cm 


■ 

10 J 


PILOT 14 




•*r 

I 


^ r 
1 


■' 


104 


STRFCV14 


" 


\r 

I 


\r 

1 


" 


103 


STFCV 14 




Y 


Y 


" 



I {. heat transfer DATA DISK -2 



No. 


FILE NAME 


DATE INJ. 


VORTEX 


DESCRIPTION 


007 


STAV2 


N 


N 


10 m/s, Stav baseline 


008 


HDAT2 


N 


N 


local h baseline. 


! 06 


TO AFC 16 


80487.1412 Y 


N 


10 m/s, M - 0.98 


107 


I D AFC 1 6 


Y 


N 


no vortex baseline 


108 


FILDT16 


Y 


N 


for 9 inj. holes F.C. 


109 


STRFC16 


Y 


N 


St/Sto 


l 10 


ST AFC 1 6 


Y 


N 


Stf(I) 


1 1 1 


STAVFC16 


Y 


N 


averaged St/Sto 


112 


TDAFC18 


80487.1452 Y 


Y 


10 m/s, M= 0.98 


113 


IDAFC18 


Y 


Y 


Vortex («> 4.79 cm 


1 14 


F1LDT18 


Y 


Y 




115 


STRFCV 1 8 


Y 


Y 


St/Sto 


1 16 


STFCV 1 8 


Y 


Y 


averaged St/Sto 
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1 17 


TDAFC20 


80487.1529 


V 

i 


Y 


10 m/s. M=0. 98- 


118 


IDAFC20 


n 


Y 


Y 


Vortex @3-52 cm 


1 19 


FILDT20 


*' 


Y 


V 

i 




120 


STRFC20 


1 1 


Y 


Y 


" 


121 


STFCV20 




Y 


Y 




122 


TDAFC22 


80487.1553 


Y 


Y 


10 m/s. M=0.98 


123 


IDAFC22 


" 


Y 


Y 


Vortex (?) 6.06 cm 


124 


^ T T TV t ^ 

r iLJi L— 


" 


i 






125 


STRFCV22 


" 


t r 

I 


V 

J. 


• 


126 


STFCV22 




v 

I 


V 

1 




127 


TDAFC24 


80987 1943 


Y 


Y 


10 m/s. M=0.98 


128 


IDAFC24 


’* 


Y 


Y 


Vortex <a 4.79 cm 


129 


FILDT24 


11 


Y 


Y 


Run for temp, profiles 


130 


STRFC24 




Y 


Y 


St/Sto 



III. TEMPERATURE PROFILES DATA DISK *1 



No. 


FILE NAME 


DESCRIPTION 


131 

132 
133 


TPR01 

TPR02 

TPR03 


10 m/s, M =0.98, heated plate, X= 1.867 m 
Vortex <» 4.79 cm , X = 2.172 m 

, X = 2.48 m 


134 

135 

136 
137 


TPR04 

TPR05 

TPR06 

TPR07 


10 m/s, M = 0.98. unheated plate, X = 1.867 m 
Vortex 4 79 cm, ,X = 2.172 m 

”, ” , X = 2.48 m 

, " , X = 1.48 m 


138 


TPR08 


10 m/s, M=0.98, heated plate, X= 1.48 m, Vt. <a 4.79 cm 


139 

140 


TPR09 

TPR010 


lOm/s, M=.98, unheated plate, X=1.48 tn. V <?> 3-52 cm 

, V (?> 6.06 cm 


141 


TPROl 1 


10 m/s, M=1.26, unheated pit., X = 1 . 4 8 m. V <a 4.79 cm 


142 


TPRO 1 2 


10 m/s, M= 0.98, unheated plate, No vortex, X=l.48 m 
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IV. TEMPERATURE PROFILES DATA DISK 3 2 



M3 

144 


TPR013 
TPRO 1 4 


10 m/s, boundary layer with fiim cooling. 0.98 

10 m/s, boundary layer heated piate only. 


i *t) 


PPROF 


10 m/s. .\rl-1.2o. v i<‘< *i./ 9 cm. Pressure cr^oe pr^me 


i 46 


CAL 


data, at X= 1.48 m 


147 


vri pivn 


velocity profiles 


1 40 


PL0TV0R 


vorticity profiles 
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